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INTRODUCTION 


The  scope  of  the  project  includes  a  number  of  separate  specialized  investigations  directed 
at  solving  problems  associated  with  yield  estimation  utilizing  short  and  long  period  body  waves. 
We  here  present  the  results  of  our  research  on  three  of  them.  In  teims  of  short  period  body  waves, 
these  problems  involve  the  effects  of  anelastic  attenuation,  pP  and  nonlinear  effects  such  as  spall 
on  seismic  amplitude  and  on  mj,.  In  a  portion  of  the  work,  we  have  developed  an  entirely  new 
approach  to  yield  estimation.  In  another  we  have  investigated  the  transportabilits  of  yield  scaling 
laws  from  site  to  site,  considering  both  U.S.  and  Soviet  testing  areas.  We  also  have  performed 
studies  on  the  effects  of  regional  and  teleseismic  Q  on  attempts  at  estimating  explosion  source 
strength. 

The  first  study  which  is  reported  on  is  an  attempt  to  generalize  the  relative  waveform  analysis 
procedure  called  intercorrelation,  to  regional  teleseismic  waves.  This  approach  has  been  used 
almost  exclusively  to  analyze  teleseismic  data  in  the  past.  There  is  a  current  need  to  characterize 
the  behavior  and  uncertainty  in  magnitude-yield  relations  down  to  body  wave  magnitudes  which 
are  as  small  as  possible.  Of  course,  the  seismic  data  with  good  signal  to  noise  ratios  for  such  small 
events  will  be  restricted  to  regional  distances.  It  has  recently  come  to  light  that  the  waveform 
in  the  western  U.S.  is  highly  stable  and  that  it  shows  a  strong  variation  with  event  size.  Intercorrelation 
has  been  used  to  measure  source  function  scaling  for  well  coupled  events  below  the  water  table 
dow  n  to  magnitude  3.9.  It  has  been  established  in  a  number  of  separate  studies  by  many  investigators 
that  the  behavior  of  pP  from  nuclear  tests  is  not  consistent  with  clastic  reflection  from  a  simple 
point  compressional  source  Direct  observation  of  isolated  pP^  appears  to  confirm  this.  Because 
of  the  relatively  low  attenuation  of  Pn  with  respect  to  teleseismic  P,  it  is  possible  to  is  late  the 
direct  P„  arrival  from  pPp.  The  parameters  of  the  source  time  function  can  be  studied  independently 
from  this  data  using  the  intercorrelation  approach.  It  was  found  that  the  relative  yield  of  events 
could  be  reliably  determined  from  the  direct  P^  wave  shapes  alone  using  thr  method. 
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The  second  investigation  we  have  undertaken  is  a  study  of  the  broad-band  nature  of  seismic 
signals  from  explosions  We  examine  both  long  period  and  short  period  typ- s  of  information 
an  attempt  to  develop  yield  scaling  laws  which  are  highly  transportable.  Along  with  the  usual 
short  period  body  waves  on  which  standard  mt-yield  curves  are  based,  we  consider  such  information 
as  long  period  P,  long  period  pS  and  near  field  signals.  We  consider  ai!  the  U.S.  test  sites  including 
those  of  GASBUGGY,  RULISON  and  FAULTLESS,  which  are  essentially  single  event  sites  and 
the  NTS  and  Amchitka  testing  areas.  We  correlate  our  results  with  those  from  the  Soviet  test  site 
at  Shagan  River. 

Though  the  new  internal  seismic  data  from  the  Soviet  Union  is  exciting  and  very  significant 
to  the  yield  estimation  process,  it  is  still  relatively  sparse  and  the  familiar  problems  of  trade-offs 
between  Q  and  explosion  source  strength  emerge.  We  have  undertaken  two  unusual  studies  of  Q, 
one  involving  regional  data  and  the  other  far  regional  to  teleseismic  data.  In  the  former,  we 
examine  Pnl  data  from  the  Soviet  station  at  GARM  and  compare  it  to  similar  data  from  the  U.S. 
station  at  Harvard.  In  the  latter,  we  consider  broad-band  data  from  Soviet  stations  KIV,  OBN 
and  ARU. 
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INTERCORRELATION  OF  REGIONAL  P„ 

The  first  study  we  will  discuss  is  an  attempt  to  take  advantage  of  some  recent  observations 
made  by  Burdick  et  al.  (1991c)  regarding  the  waveform  of  the  onset  of  regional  It  appears 
clear  that  these  waveforms  are  very  directly  comparable  to  those  of  teleseismic  P.  Rather  than 
being  true  head  waves,  these  arrivals  seem  to  be  turning  rays  in  the  uppermost  mantle  above  the 
low  velocity  zone.  At  short  times,  they  can  be  represented  as  a  simple  convolution  of  a  source,  Q 
and  instrument  operator  just  as  teleseismic  short  period  P  can.  However,  because  they  never 
penetrate  the  low  velocity  zone,  the  effect  of  attenuation  is  tremendously  reduced.  Instead  of  t* 
having  a  value  of  approximately  1.0  s,  as  it  does  for  teleseismic  P  from  NTS,  it  has  a  a  value  close 
to  0.1s;  an  order  of  magnitude  change  in  the  exponential.  This  allows  the  direct  P  to  be  separated 
from  the  pP  through  simple  time  domain  windowing. 

Before  discussing  the  regional  P^  yield  scaling  behavior,  we  will  review  the  procedure  known 
as  waveform  intercorrelation.  It  is  a  method  which  allows  the  reduction  of  the  observations  of 
many  events  at  a  common  test  site  to  source  models  for  the  events.  The  source  models  are 
parameterized  mathematical  functions  which  describe  the  reduced  displacement  potential  and,  for 
teleseismic  data,  a  representation  for  the  free  surface  interaction  and  perhaps  even  spall.  The  RDP 
representation  is  typically  an  explicit  function  of  yield  with  the  free  parameters  defined  in  terms 
of  >  ield  scaling  laws.  The  changes  in  the  shape  of  the  RDP  for  events  at  a  common  site  are  thus 
related  to  changes  in  yield  The  key  advancement  in  this  study  is  that  the  free  surface  parameters 
can  be  neglected.  The  most  important  goal  is  to  attempt  to  estimate  yield  directly  from  the  evolution 
of  the  shape  of  the  P^  onset. 

The  classic  intercorrelation  procedure  is  illustrated  in  Figure  1.  Each  seismogram  from  a 
given  event  is  convolved  with  the  source  function  and  a  spike  ttain  representint’  P  and  pP  from 
another  event.  The  seismograms  from  the  other  event  are  convolved  with  the  P+pP  spike  train 
and  parameterized  RDP  from  the  first  event.  The  two  intercorrelation  waveforms  at  each  station 
are  then  analytically  compared  to  each  other  and  the  source  functions  that  produce  the  best  w  aveform 
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Figure  1.  An  Illustradon  of  (he  standard  intercorrelation  procedure.  The  left  column  shows  two 
teleseismic  short  perHvi  P  waveforms  from  a  common  test  site  at  a  common  station.  The  center 
column  shows  source  estimates  for  each  event  paired  with  the  signals  from  the  opposite  event.  The 
source  estimates  are  described  in  terms  of  both  a  source  time  function  and  pP  reflection.  The 
cross-paired  sources  and  records  are  convolved  and  the  resultants  equlai  red  in  a  crid  search  operation 
using  an  analytic  norm  The  results  foi  this  case  are  show  n  on  the  right  In  the  new  prcxredure  used 
here,  the  pP  phase  is  v.  indowcd  out  and  only  the  source  parameters  arc  searched  for. 
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match  are  found.  The  procedure  is  applied  to  all  possible  matching  stations  simultaneously.  In 
th  >  instance,  these  are  tv.o  Novaya  Zemlya  events  as  observed  at  station  ATL.  The  equalization 
of  the  two  dissimilar  waveforms  on  the  left  after  the  intercorrelation  procedure  is  very  good  as 
shown  on  the  right. 

To  incorporate  differences  in  event  yields,  it  is  necessary  to  adopt  a  parameterized  description 
of  the  explosion  source.  Such  parameterizations  are  abundant  in  the  literature  (Murphy,  1977; 
von  Seggern  and  Blandford,  1972;  Helmberger  and  Hadley,  1981).  Here  we  adopt  the  Murphy 
(1977)  source  representation  and  yield  scaling  though  we  emphasize  that  the  other  two  are  equa.'y 
valid  and  would  perform  as  well  if  need  be.  The  intercorrelation  waveforms  for  each  station  pair 
are  typically  compared  using  one  or  two  norms.  The  waveform  norm  is  given  by 

-  -1(1  -CCC) 
n 

CCC  is  the  normalized  cross  correlation  coefficient  and  n  is  the  number  of  stations.  This  norm  is 
most  sensitive  to  differences  in  zero  crossing  times  and  is  often  used  in  waveform  inversion  studies. 
Absolute  amplitude  information  is  typically  retained  using  a  standard  least  squares  norm  if  need 
be,  but  in  this  investigation  we  have  utilized  only  the  waveform  norm. 

Over  the  past  few  decades  the  state  of  digital-broad-band  seismic  recording  has  gone  through 
a  major  revolution.  Fortunately  for  the  endeavor  of  test  ban  treaty  monitoring,  many  of  the  highest 
quality  stations  were  located  in  the  western  U.S.  where  they  recorded  many  NTS  nuclear  events. 
Figure  2  is  a  map  of  the  stations  from  which  we  have  assembled  the  data  base  we  have  used  in 
th  is  work.  The  stations  shown  were  installed  over  a  substantial  period  of  time  and  illustrate  the 
development  of  digital  seismic  recording.  The  ones  which  were  installed  first  were  DWWSSN 
stations  ALQ  and  JAS.  Recording  was  and  is  made  in  several  pass-bands  and  is  available  on  the 
standard  Network  Day  Tapes.  The  four  LLNL  stations,  MNV,  FLK,  KNB  and  LAC  represent  a 
major  upgrade  in  recording  capability  because  they  are  broad-band  with  a  wide  dynamic  range. 
(We  note  that  there  were  important  analog  stations  at  these  same  sites  decades  in  the  past.)  The 
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Figure  2.  The  western  U.S.  digital  array  as  defined  Tor  the  purpose  of  this  study.  MNV,  KNB, 
LAC  and  ELK  are  LLNL  stations,  .ALQ  amd  JAS  are  DWWSSN  st.ilions,  and  PAS  and  PFO  are 
university  run  stations.  The  open  circles  are  the  location  of  earthquakes  in  tlie  discrimination  data 
base. 
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two  stations  PAS  and  PFO  have  state-of-the-art  Streckeisen  instruments  which  are  the  same  as 
th(  se  used  in  the  IRIS  net.  The  important  point  in  terms  of  his  study  is  that  they  all  record  regional 
P„  and  have  provided  an  ample  data  base  for  this  investigation. 

The  concept  of  using  the  intercorrelation  methodology  to  measure  yield  from  the  onset  of 
regional  Pn  is  not  .vithout  foundation.  In  a  recent  work,  Burdick  et  al.  (1991c)  showed  that  this 
onset  was  simple,  stable  and  easy  to  interpret  in  terms  of  the  standard  concepts  used  to  analyze 
teleseismic  short  period  P  waves.  Figure  3  shows  some  of  the  typical  data  from  the  Burdick  et  al. 
(1991c)  study.  Observed  deconvolved  velocity  waveforms  from  4  Yucca  Flat  events  as  recorded  at 
digital  station  MNV  are  shown  in  solid  line.  Almost  1.5  orders  of  magnitude  are  represented 
between  the  largest  and  smallest  events.  The  frequency  difference  between  them  is  clear  even  in 
the  time  domain.  Synthetics  computed  using  the  Murphy  (1977)  source  and  an  assumption  of 
elastic  pP  behavior  are  shown  as  dashed  lines.  The  observed  and  predicted  pP  arrivals  are  indicated 
by  arrows.  The  observed  pPn  arrivals  are  late  and  close  to  the  elastic  predictions  in  size,  but  in 
any  case,  in  this  representation  it  is  simple  to  window  out  the  distinct  direct  Pn  arrival.  In  passing, 
we  note  that  these  observations  lend  additional  credence  to  the  intercorrelation  results  as  opposed 
to  spectral  averaging  results.  We  also  note  that  the  Murphy  (1977)  predictions  of  change  in  frequency 
content  between  the  largest  and  smallest  event  are  not  satisfactory.  The  predicted  changes  in  time 
function  are  less  extreme  than  the  observed.  The  data  available  from  LLNL  is  limited,  but  in 
resiewing  it  we  did  find  signals  from  eleven  Yucca  Flat  events  below  the  water  table.  Figure  4 
compares  the  waveforms  at  common  LLNL  stations  for  the  smallest  and  the  largest,  BORREGO 
(m^,  =  3  9)  and  JORNADA  (mj,  =  5.9).  The  change  in  frequency  content  at  all  stations  is  dramatic. 
Under  each  complete  Pn  trace  is  a  windowed  trace  where  the  window  is  a  trapezoid  with  an  0.5  s 
lead  a  1.0  s  level  beginning  at  estimated  onset  time  and  an  0..*'  second  fall.  There  are  clearly 
waveform  differences  between  the  stations,  particularly  at  ELK,  but  the  principle  of  intercorrelation 
is  to  characterize  the  changes  in  waveshape  at  fixed  stations  allowing  for  the  possibility  of  complex 
path  effects.  The  window  should  allow  good  spectral  measurements  up  to  frequencies  of  I  hz. 
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Pn  from  YUCCA  of  MNV 

(Ground  Velocity  for  Turning  Roy) 


Figure  3.  An  illustration  of  the  Pn  waveform  observations  which  motivate  this  study.  Observed 
waveforms  are  shown  as  solid  lines  and  synthetics  as  dashed.  The  pPn  is  a  downward  pulse  as 
indicated  by  an  arrow  in  both  the  data  and  the  synthetics.  The  pPn  is  delayed  in  the  observations 
and  can  be  windowed  out.  This  permits  intercorrelation  analysis  of  the  isolated  direct  arrival. 
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f  igure  4.  An  illustration  of  the  windowing  procedure.  The  top  trace  shows  the  initial  pulse.  It  is 
windowed  with  a  trapa7()id  with  an  0.5  s  rise,  an  0.75  s  level  time  commencing  at  estimated  arrival 
time  and  an  0.5  s  fall.  1  he  removal  or  at  least  strong  reduction  of  the  pPn  is  apparent. 
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The  procedure  for  intercorrelation  has  been  presented  in  many  previous  reports  (Burger  et 
a!.,  1986,  Lay  et  al.,  1984a,  Lay,  1985,  1991).  As  discussed  above,  the  basic  idea  is  to  choose 
reference  event  and  an  estimated  source  function  for  that  event.  If  the  signal  includes  pP,  then 
pP  parameters  must  also  be  included,  but  in  this  case,  only  direct  Pn  is  present.  We  selected 
BORREGO  and  used  the  Murphy  (1977)  source.  The  procedure  is  then  to  find  a  source  function 
for  each  other  event  which  minimizes  the  differences  between  waveforms  when  a  cross-convolution 
is  performed.  Since  the  depth  of  these  events  is  known,  the  only  free  parameter  is  yield,  and  we 
simply  stepped  through  yield  values  to  find  those  for  the  non-reference  events  which  optimized 
the  waveform  norm.  The  yield  of  BORREGO  was  estimated  to  be  0.7  kt  from  the  scaling  relation 
of  Yacoub  (1984)  based  on  WWSSN  records  (Network  A  A)  of  Pahute  and  Yucca  events  below  the 
water  table.  Figure  5  shows  an  additional  unprocessed  example  for  BORREGO  and  BOUSCHET 
which  are  also  Yucca  Flat  events  below-  the  water  table.  This  pairing  of  events  is  quite  remarkable 
since  they  had  identical  depths.  The  shift  in  frequency  content  is  quite  clear,  and  because  the 
depths  are  identical,  it  is  almost  certainly  due  to  yield  scaling.  Figure  6  shows  the  results  of  the 
intercorrelation  of  those  two  events.  The  fits  are  quite  remarkable  in  this  example.  It  should  also 
be  noted  that  we  are  cross  correlating  waveforms  from  events  differing  by  two  orders  of  magnitude 
in  size.  The  match  of  the  complexity  at  ELKO  is  worth  special  recognition,  and  it  would  seem 
to  indicate  that  the  intercorrelation  approach  accounts  for  complex  site  effects  as  it  was  designed 
to  do. 

Figures  7  and  8  show  similar  initial  and  intercorrelated  results  for  Yucca  events  above  the 
water  table.  The  reference  event  was  selected  as  CORREO  and  the  test  event  was  TENAJA.  In 
general,  it  seems  best  to  select  the  smallest  events  for  reference,  perhaps  because  they  have  the 
richest  frequency  content.  That  is  to  say  that  the  for  the  largest  events,  the  high  frequency 
information  is  washed  out.  The  complexity  of  the  ELKO  station  is  still  apparent  and  the  equalized 
records  for  that  station  are  comparable  for  events  above  and  below  the  water  table.  The  assembly 
of  a  data  base  for  Pahute  Mesa  events  below  the  water  table  is  problematical  since  onlv  a  few  have 
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Figure  5.  These  are  sample  observations  from  one  of  the  events  in  the  data  base  of  Yucca  Flat 
events  below  the  water  ttible.  The  reference  event  on  the  right  (BORREGO)  had  a  depth  of  564  m 
but  a  magnitude  of  only  3.9.  BOUDCFIET  had  a  magnitude  of  5.7  at  the  same  depth.  The  shift  in 
frequency  content  due  to  yield  dependent  variation  in  the  source  function  is  clear.  Since  it  can  be 
readily  observed  it  can  be  quantified  through  intercorrelation  analysis. 
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SAMPLE  INTERCORRELATION 

BOUSCHET  BORREGO 

Yucco 

Below  W.T. 


Figure  6.  These  are  the  intercorrelated  signals  shown  in  unprocessed  form  in  Figure  5.  The 
waveforms  have  been  clearly  equalized.  The  yield  of  the  reference  event  was  estimated  at  0.7  kt 
from  a  standard  NTS  mb  -  yield  relationship.  The  yield  of  BOUSCHET  was  estimated  at  80  kt. 
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F  igure  7.  These  are  illustrative  waveforms  from  the  data  base  of  Yucca  Flat  events  above  the  water 
table.  1  he  events  are  shallower  than  those  in  the  previous  example,  .so  windowing  out  of  pP  is  less 
effective.  The  reference  event  selected  for  the  analysis  of  this  data  base  was  CORREO  with  a  depth 
of  335  m  and  a  magnitude  of  4.8.  In  this  example  TENAJA  has  a  depth  of  357  m  and  a  magnitude 
of  4.5. 
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Figure  8.  These  are  the  intercorrelated  signals  shown  in  unprocessed  form  in  Figure  7.  The 
waveforms  have  again  been  effectively  equalized  using  only  the  source  function  in  the  intercor¬ 
relation.  The  yiels  of  the  reference  event  was  estimated  at  9.4  kt.  The  yield  of  TENAJA  was 
estimated  at  10  kt. 
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been  detonated  since  the  threshold  test  ban  agreement  of  ISO  kt.  In  the  data  base  available  from 
th(  LLNL  digital  stations,  only  two  such  events  were  recorded.  Id  perform  a  partially  legitimate 
test  of  our  approach,  we  again  selected  BORREGO  (Yucca  Flat  below  the  water  table)  as  a  reference 
event.  Figure  9  shows  the  raw  signals  and  Figure  10  shows  the  intercorrelation  for  event  TOVVANDA. 
The  equalization  of  the  waveforms  is  remarkable  and  the  unusual  complexity  of  the  ELKO  record 
persists.  Finally,  Figures  1 1  and  12  show  samples  of  the  raw  data  and  intercorelations  for  Pahute 
events  above  the  water  table.  The  reference  is  KAPPELI  and  the  test  event  is  HOSTA.  The 
procedure  again  appears  successful. 

As  discussed  above,  we  applied  intercorrelation  analysis  to  all  the  events  available  in  our 
data  base  to  predict  the  yield  scaling  behavior  and  to  compare  it  to  expected  results.  We  review 
here  exactiv  how  the  analysis  proceeds  once  the  relative  yield  has  been  determined.  There  are  two 
basic  goals  in  the  structure  of  the  analysis.  The  first  is  simply  to  estimate  a  yield  for  the  master 
event.  The  second  and  much  more  important  goal  is  to  loosely  relate  our  results  to  the  widely 
accepted  rules  for  yield  scaling.  There  is  little  debate  regarding  the  mt-yield  scaling  relationships 
for  events  larger  than  100  kt  at  NTS.  From  several  possible  choices,  we  selected  the  scaling  law 
of  ^  acoub  (1984).  We  then  accept  either  the  ISC  or  LLNL  body  wave  magnitude  (ISC  preferred) 
for  the  master  event.  We  then  derive  a  yield  for  the  master  event  from  the  Yacoub  (1984)  law. 
From  this  point  forward,  the  procedure  is  completely  automated  and  depends  only  on  the 
intercorrelation  analysis.  Figure  13  shows  the  results  for  Yucca  flat  events  below  the  water  table. 
The  selected  master  event  was  BORREGO.  The  solid  line  in  the  graph  is  the  Yacoub  (1984)  curve 
and  thus  the  master  event  falls  directly  on  it.  The  other  events  scatter  evenly  about  the  predicted 
curve.  One  of  the  events  is  JORNADA  which  is  compared  w  th  BORREGO  in  Figure  4.  The 
implication  is  clear.  The  strong  change  in  frequency  content  of  direct  P^  in  that  figure  is  directly 
related  to  yield.  Furthermore,  the  large  event  scaling  law  appears  to  hold  down  to  events  as  small 
as  BORREGO. 
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Figure  9.  The  waveforms  on  the  left  come  form  an  event  below  the  water  table  at  Pahute  Mesa. 
Only  t^»'0  events  of  almost  identical  size  were  available  in  the  data  base  so  BORREGO  (Yucca 
below  water  table)  was  used  as  a  reference  as  in  Figure  5. 
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SAMPLE  INTERCORRELATION 

TOWANDA  BORREGO 

Pahute 

Below  W.T. 


ELK 


Figure  10.  These  are  the  intercorrelated  signals  shown  in  unprocessed  form  in  Figure  9.  The 
waveforms  have  again  been  effectively  equalized  using  only  the  source  function  in  the  intercor- 
rclation.  The  yield  of  the  reference  event  was  estimated  at  0.7  kt  and  of  TOWANDA  at  24  kt. 
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Figure  1 1.  These  waveforms  are  from  Pahute  events  above  the  water  table.  The  magnitude  of 
reference  event  KAPPELI  is  5.2  and  of  test  event  HOSTA  is  5.6. 
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Figure  12.  These  are  the  intercorrelated  signals  shown  in  unprocessed  form  in  Figure  10.  The 
waveforms  have  again  been  effectively  eqqualized  using  only  the  source  function  in  the  intercor¬ 
relation.  The  yields  of  the  reference  event  was  estimated  at  35  kt  and  of  HOSTA  at  24  kt. 
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Ligiirc  1  3.  A  comparison  of  yield  sc;iling  predictions  based  on  only  intercorrelation  analysis  of  Pn 
compared  to  a  standard  scaling  law  for  NTS.  In  this  case,  the  data  base  included  Yucca  events 
below  the  water  table.  The  reference  event  was  BORREGO. 


Figure  14  shows  the  results  for  Yucca  Flat  events  above  the  water  table.  In  this  case,  the 
ref>  rence  event  was  chosen  as  COR  REO  which  is  approximately  in  the  middle  of  the  m^  distribution. 
The  reason  for  this  choice  is  that  CORREO  has  the  most  complete  recorded  data  base,  though  it 
is  true  that  the  smallest  master  e\ents  tend  to  give  the  best  intercorrelation  results.  CORREO  is 
nea  the  top  of  the  population  as  it  distributes  about  the  standard  yield-scaling  line.  The  remarkable 
fact  is  that  the  slope  of  the  population  follows  very  closely  the  prediction  of  the  Yacoub  (1984) 
line.  There  is  an  open  question  as  to  whether  there  are  changes  in  slope  of  this  line  as  events 
become  very  small.  The  evidence  here  is  that  the  slope  does  not  change  for  events  with  yields  as 
low  as  about  2  kt.  Also  note  that  the  distribution  of  events  is  much  more  even  than  in  Figure  13. 
There  seems  to  be  little  difference  in  behavior  for  events  above  and  below  the  water  table. 

Figure  15  shows  the  results  for  Pahute  events  below  the  water  table.  As  noted  previously, 
there  have  been  few  such  tests  in  recent  history.  BORREGO  from  beneath  the  water  table  was 
arbitrarily  selected  as  the  master.  The  Yacoub  (1984)  prediction  holds  reasonably  well  although  a 
small  increase  in  slope  may  be  indicated  by  the  data  points.  Figure  16  shows  the  results  for  Pahute 
Mesa  events  above  the  water  table.  The  master  event  was  chosen  .as  KAPPELI.  In  this  instance, 
there  appears  to  be  an  indication  of  a  decrease  in  slope  in  the  data  points.  However,  the  distribution 
in  magnitudes  is  quite  limited.  Clearly,  the  analysis  of  more  data  with  a  greater  range  of  magnitudes 
would  be  desirable  in  the  future. 

Discu.sslQn:  The  importance  of  these  discoveries  in  terms  of  modern  treaty  monitoring  is  very 
substantial.  Though  there  is  clearly  the  need  for  much  more  investigation,  the  onset  method 
of  estimating  yields  represents  a  fourth  new  seismic  approach  for  doing  so,  and  one  that  can  be 
used  for  very  small  events.  The  three  classic  methods  are  to  examine  m^-yield  scaling,  Mg-yield 
scaling  or  Mg-yield  scaling.  The  difference  with  the  P„  onset  approach  is  that  it  relies  on  wave 
shape  changes  as  a  function  of  yield  as  opposed  to  amplitude  information,  though  the  latter  could 
certainly  be  incorporated  into  the  procedure  in  the  future.  Ms  scaling  is  generally  simply  a  measure 
of  the  dependence  of  20  second  spectral  amplitude  on  yield.  Complete  waveform  information  is 
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Figure  14.  A  comparison  of  yield  scaling  predictions  based  on  only  intercorrclation  analysis  of  Pn 
compared  to  a  standard  scaling  law  for  NTS.  In  this  case,  the  data  base  included  Yucca  events 
above  the  water  table.  The  reference  event  was  CORREO. 
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Figure  13.  A  comparison  of  yield  scaling  predictions  based  on  only  intercorrelation  analysis  of  Pn 
compared  to  a  standard  scaling  law  for  NTS.  There  were  only  two  events  available  from  Pahute 
below  the  water  table.  The  reference  event  used  was  BORREGO  (Yucca  below  the  water  table). 
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Figure  16.  A  comparison  of  yield  scaling  predictions  based  on  only  intercorrelation  analysis  of  Pn 
compared  to  a  standard  scaling  law  for  NTS.  In  this  case,  the  data  base  included  Yucca  events 
above  the  water  table.  The  reference  event  was  KAPPELI.  In  this  instance  there  appears  to  be 
some  difference  in  slope  between  the  intercorrelation  measurements  and  the  standard  law. 
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not  utilized.  M®  or  moment  scaling  does  utilize  much  more  of  the  long  period  information,  but 
bot  of  these  techniques  suffer  from  serious  difficulties  at  magnitudes  less  than  4.0.  Body  wave 
magnitude  does  depend  on  a  single  period  measurement  in  that  it  scales  as  Log(A/T),  but  it  hardly 
utilizes  complete  waveform  information.  Furthermore,  m^-yield  scaling  is  the  most  widely  accepted 
approach  to  yield  estimation,  and  in  those  magnitude  ranges  where  both  m^  values  and  waveforms 
are  available  both  methods  give  the  same  results. 

There  are  many  directions  which  need  to  pursued  in  the  future  regarding  yield  estimation 
from  P„  onset  waveforms.  The  generality  of  the  approach  at  alternate  test  sites  remains  a  significant 
issue,  though  the  physics  of  the  approach  are  very  simple  and  there  is  every  reason  to  believe  it 
will  be  transportable.  The  variability  of  the  results  with  respect  to  events  above  and  below  the 
water  table  at  Pahute  Mesa  is  an  important  issue,  but  it  can  be  investigated  further  by  considering 
additional  NTS  events. 
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BROADBAND  SOURCE  MODELS  FOR  U.S.  UNDERGROUND  NUCLEAR  EXPLOSIONS 
Intrnrfuctinn:  Recent  broadband  studies  of  earthquakes  indicate  considerable  complexity  and 
non-uniformity  in  source  characteristics.  Modern  source  descriptions  are  expressed  in  terms  of 
seismic  moment  and  asperity  distribution  on  the  fault  surface.  The  latter  is  best  established  by 
studies  of  local  strong-motion  observations  while  the  former  can  be  obtained  from  teleseismic 
modeling  of  long-period  body  waves  and  surface  waves.  The  recent  deployment  of  the  new  IRIS 
systems  consisting  of  the  Weilandt-Strekeisen  sensors  and  24  bit  Quanterra  loggers  allows  the  entire 
frequency  band  to  be  recorded  and  modeled  locally  (e.g.,  Dreger  and  Helmberger,  1990).  Since 
smaller  events  can  not  be  seen  teleseismically,  the  new  data  systems  become  essential  in  studying 
both  earthquake  and  explosion  sources  and  in  establishing  techniques  for  discrimination.  The 
broadband  nature  of  underground  nuclear  sources  is  the  key  question  in  discrimination.  Furthermore, 
the  detailed  variations  in  these  source  descriptions  for  various  test  sites  in  the  U.S.  and  in  foreign 
environments  are  also  important.  We  will  address  these  issues  with  respect  to  U.S.  explosions  and 
attempt  to  establish  some  useful  kinemetric  source  descriptions. 

Source  models  and  scaling  relationships  for  underground  explosions  have  been  studied  for 
several  decades,  but  generally  with  an  emphasis  on  comparing  teleseismic  data  sets  consisting  of 
short-period  wt’s  and  long-period  surface  wave  A/s's.  Classical  amplitude  measures  on  the  F-wave 
such  as  can  only  take  account  of  complicating  factors  such  as  attenuation,  pP,  and  initial  source 
histories  in  a  crude  fashion.  Similar  measures  on  surface  waves  (Afs)  suffer  from  their  own 
complicating  factors  such  as  effective  attenuation,  tectonic  contamination  and  source  coupling.  It 
has  been  apparent  for  many  years  that  these  two  yield  estimators  show  distinct  regional  behavior, 
presumably  caused  by  differences  in  the  above  factors.  Figure  17a  displays  measurements  taken 
from  a  sample  of  NTS  and  SRS  events.  The  regression  lines  for  these  two  regions  indicate  distinctly 
different  trends.  The  separation  of  these  populations  is  well  established  at  the  larger  yields,  as 
demonstrated  in  many  studies  (Sykes  and  Cifuentes,  1984;  Given  and  Mellman,  1986). 
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Figure  17.  (a)  Plots  of  vs.  Mg  regression  curves  for  Shagan  River,  USSR,  and  NTS  explosions 
(after  Marshall  el  al.,  1979).  (b)  Same  regression  curves  as  in  (a)  shown  with  Soviet  PNE’s  and 
off-NTS  explosions  CANNIKIN,  FAULTLESS,  GASBUGGY,  LONGSHOT,  MILROW,  and 
RULISON  (indicated  by  the  first  letters  of  their  names.) 
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Part  of  the  scatter  displayed  in  Figure  17a  is  caused  by  random  effects  and  part  is  caused 
by  deterministic  factors  which  can  be  modeled  and  are  correctable.  Separating  these  two  el  iects 
is  fundamental  to  higher  resolution  attempts.  Figure  17b  shows  the  results  of  wt  vs.  A/g  comparisons 
for  a  sample  of  peaceful  nuclear  explosions.  They  show  a  great  deal  of  variation  which  we  can 
address  to  help  isolate  the  deterministic  elements  of  the  source  process  and  establish  source  properties. 

Numerous  papers  have  been  written  on  why  wt  and  Wg  values  f'om  Amchitka  shots  are 
different  from  those  •'T  Pahute.  One  of  the  earliest  reports  was  given  by  Von  Seggern  (1972)  who 
investigated  the  differences  between  BOXCAR  (Pahute)  and  MILROW  (Amchitka),  two  events  of 
nearly  the  same  yield.  The  for  BOXCAR  is  nearly  0.3  units  less  than  that  of  MILROW  while 
the  A/g  for  BOXCAR  is  approximately  0.5  units  larger  than  that  of  MILROW.  He  concluded  that, 
since  teleseismic  P-waves  from  NTS  are  lower  frequency  than  those  from  Amchitka,  greater 
attenuation  under  NTS  was  the  reason  for  the  differences  in  m^.  He  gave  no  explanation  for  the 
differences  in  A/g  between  the  two  events.  Similarly,  the  recent  Joint  Verification  Experiment, 
involving  shots  of  the  same  yield  in  the  U.S.  and  USSR,  showed  a  similar  mb  offset.  This  is 
generally  believed  to  be  caused  by  a  difference  in  the  attenuation  beneath  the  two  test  sites. 
However,  the  Afg’s  from  Shagan  River  events  which  are  the  least  contaminated  by  tectonic  release 
(as  indicated  by  the  absence  of  Love  waves)  are  roughly  0.3  to  0.5  units  less  than  NTS  events  for 
events  near  the  I50kt  testing  limit  (Stevens,  1986).  Sykes  and  Cifuentes  (1984)  used  these  same 
events  to  argue  for  Soviet  compliance.  These  discrepancies  suggest  that  the  long-period  source 
levels  sampled  by  .Vfg  relative  to  the  short-period  levels  sampled  by  mb  are  site-dependent  and, 
thus,  these  differences  can  be  used  to  establish  some  working  broadband  source  models. 

The  approach  followed  here  is  similar  to  that  of  Lay  ei  al.  (1984b).  We  assume  a  convenient 
modified  Haskell  source  representation  given  by 

V(/)-  V.^1  1  ♦  A/*  ^  -  S(A:t)'|j 
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where  the  RDP  is  defined  by  V(t)/^.  The  reduced  time  •tisT-(-(^/a)  where  R  is  the  distance 

b<  tween  the  source  and  receiver  and  a  is  the  velocity.  The  study  by  Mueller  and  Murphy  (1971) 
established  the  basic  scaling  laws  relating  the  constants  K  and  to  yield.  Their  formalism  is 
easily  adapted  to  the  above  RDP.  The  parameter  K  is  directly  related  to  the  corner  frequency. 
The  parameter  B  controls  the  overshoot  which  provides  the  means  to  uncouple  the  short-period 
signals,  used  for  the  measurements,  from  the  static  level,  which  controls  the  Mg  measurement 
(Figure  18).  Thus,  allowing  B  to  vary  provides  the  extra  freedom  in  modeling  extended  data  sets. 
This  study  will  address  the  determinations  of  the  three  constants,  K,  B  and  as  a  function 

of  energy  (yield)  for  the  various  U.S.  test  sites  including  the  PNE’s.  We  will  begin  with  the  simplest 
data  set,  Amchitka,  followed  by  the  most  complex,  NTS. 

Scaling  Relalionshins  Amchitka:  The  data  set  for  Amchitka  consists  of  near-in  strong  motion 
seismograms,  and  teleseismic  seismograms  containing  body  waves  and  surface  waves  (see  Figure 
19).  In  earlier  efforts  Burdick  et  al.  (1984)  and  Lay  et  al.  (1984a)  established  the  linkage  between 
these  types  of  near-in  data.  Figure  19a,  and  the  teleseismic  short-period  amplitudes.  Their  source 
models,  while  fitting  the  short-period  data  very  well,  do  not  fit  the  relative  amplitudes  of  the 
Rayleigh  waves  given  in  Figure  19b.  They  predict  CANNIKIN/MILROW  long-period  amplitude 
ratios  of  about  3  assuming  B=\.  The  observed  ratio  is  near  6.  Note  that  the  Amchitka  events  have 
negligible  Love  waves,  so  we  would  not  expect  to  see  any  tectonic  effects  on  the  Rayleigh  waves. 
The  amplitude  ratio  of  the  observed  long-period  P-waves  between  these  two  events  is  about  2.8 
which  is  not  consistent  with  the  short-period.  The  Mg’s  from  the  Amchitka  shots  are  3.9 
(LONGSHOT).  4.9  (MILROW),  and  5.7  (CANNIKIN)  (Liebermann  and  Basham,  1971;  Willis  et 
al.,  1972)  and  fit  the  well  known  equation 

Ms-Logy  +  2 

which  is  an  often  quoted  result. 
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Figure  18.  Reduced  displacement  potentials  (RDP)  for  a  fixed  K  (=6)  and  and  various  values 
of  B  (1.00,  0.75,  0.50,  0.25,  0.00).  This  figure  demonstrates  how  B  controls  the  degree  of  overshoot 
in  RDP’s. 
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Figure  19.  Shown  in  the  top  panel  is  Amchitka  Island  and  records  from  events  CANNIKIN  and 
MILROW.  Three  component  velocity  seismograms  are  presented  along  with  the  slant  distances 
from  event  to  station.  Bottom  panel  shows  long-period  vertical  Pand  Rayleigh  waves  for  CANNIKIN 
and  MILROW  at  station  DUG  (Dugway,  Utah).  For  MILROW'  the  body  wave  is  larger  but  for 
Cannikin  the  surface  wave  is  larger. 
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Lay  et  al.  (1984b)  allowed  B  to  vary  and  resolved  this  discrepancy.  We  will  review  some  of 
these  efforts  and  modify  the  final  scaling  laws  to  better  match  absolute  amplitude  levels. 
Modeling  Near-In  Data:  About  six  good  sets  of  strong  motion  records  of  the  type  displayed  in 
Figure  19a  are  available  for  both  MILROW  and  CANNIKIN.  A  crustal  model  was  derived  to  fit 
the  travel  times  of  these  records  and  many  others.  These  crustal  structures  were  further  refined 
by  matching  the  entire  recorded  waveform  with  the  synthetic  waveform  (Burdick  el  al.,  1984). 
The  most  stable  portion  of  these  observed-synthetic  matches  is  the  initial  few  seconds  of  motion 
as  displayed  in  Figure  20a.  The  synthetic  consists  of  the  diving  P  followed  by  the  opposite  in 
polarity  pulse,  pP  (see  Vidale  and  Helmberger,  1987,  for  details  and  the  effects  of  2D  structures). 
However,  since  pP  cancels  P  it  becomes  difficult  to  resolve  all  three  parameters  B,  K,  and  V. 
because  of  the  4'.  vs.  B  trade-off,  as  displayed  in  Figure  20b.  For  instance,  using  K  =  6,  a  source 
function  with  V.  -  7.3and  B  -  0.5  fits  the  data  just  as  well  as  a  source  function  with  Y.  -  4. Sand 
B  -  1.  Holding  fi=l  allows  an  estimate  of  V.  for  the  study  events.  These  VJs  can  be  compared 
with  large  teleseismic  short-period  data  sets  to  establish  a  realistic  estimate  of  i‘„  and  in  this  case, 
ti-0.9(as  reported  by  Burdick  et  al.  1984).  Adding  the  intercorrelation  procedure  (Lay  et  al., 
1984b)  allows  source  strengths  of  other  events  such  as  LONGSHOT  to  be  estimated  and  a  scaling 
law  developed  based  on  short-period  signals  (see  Figure  2Ia.)  Lay  et  al.  {1984a)  fixed  B  at  1  for 
MILROW  and  adjusted  the  B'%  of  LONGSHOT  and  CANNIKIN  to  match  the  Rayleigh  wave 
differentials  (Figure  21b).  A  further  modification  is  possible  if  we  can  establish  the  absolute 
long-period  level. 

Three  types  of  data  are  available  for  this  purpose,  namely  long-period  F-waves,  the  phase 
pS  and  the  direct  modeling  of  the  Rayleigh  waves.  The  long-period  P- waves  prove  disappointing 
because  the  interactions  of  pP  with  /’  reduces  the  sensitivity  to  B.  Thus  ratios  (SPZ/LPZ)  of 
short-period  vertical  P-wave  amplitudes  to  long-period  vertical  Rayleigh  wave  amplitudes  depend 
mostly  on  the  values  of  /*  (see  Figure  22).  The  average  observed  ratio  SPZ/LPZ  using  28  long-period 
records  is  0.65,  and  we  again  obtain  a  i  \  near  I.O.  The  phase  pS  appears  the  most  promising,  since 
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Figure  20.  A  suite  of  source  models  for  event  CANNIKIN  (top  panel).  The  H'.  determined  from 
the  first  swing  in  the  waveform  is  printed  to  the  left  of  each  synthetic.  Bottom  panel  shows  the 
trade-off  of  y.  and  5  for  CANNIKIN  for  a  specific  K. 
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Figure  21.  Reduced  displacement  potentials,  (a)  based  on  short  period  signals,  (b)  from  broadband 
data:  fixing  B-\  for  MILROW,  and  adjusting  the  B'%  of  LONGSHOT  and  CANNIKIN  to  match 
Rayleigh  wave  differentials,  (c)  from  broadband  data:  utilizing  long-period  P-waves,  the  phase 
pS,  and  Rayleigh  waves.  (Figure  modified  from  Lay  el  al.,  1984.] 
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Amplitude  Ratios  for  CANNIKIN 


Figure  22.  Amplitude  ratios  of  short  period  vertical  (SPZ)  to  long  period  vertical  (LPZ)  for  event 
CANNIKIN.  The  average  observed  ratio  using  28  LPZ  records  is  0.65,  indicated  by  the  solid  line. 
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these  waveforms  are  not  subjected  to  any  strong  interference  such  as  the  P-pP  waveforms.  Some 
examples  of  the  21  observations  available  are  given  in  Figure  23  However,  one  complication  with 
this  phase  is  the  S-P  conversion  at  the  Moho  beneath  the  receiver.  Note  the  small  positive  pulse 
preceding  the  negative  direct  pS  arrival.  There  can  also  be  problems  associated  with  later  arrivals, 
the  so-called  JF-coupled  PL  waves.  Other  complications  can  occur  at  the  source  region  caused 
by  tectonic  release  and  possible  spall  interaction.  These  problems  are  discussed  at  length  by  Cohee 
and  Lay  (1988)  in  their  study  of  pS  from  some  Novaya  Zemlya  events. 

The  synthetics  displayed  in  Figure  23  were  calculated  assuming  that  the  elastic  conversion 
F  to  5  is  -0.456  for  a  ray  parameter  of  0.12  sec/km.  This  is  consistent  with  the  value  used  in 
modeling  the  phase  pS  in  earthquake  studies  (Langston  and  Helmberger,  1975).  Note  that  the 
actual  reflection  point  occurs  outside  the  region  of  strongest  spall  (a  distance  of  roughly  1.5  km 
for  CANNIKIN)  because  of  the  change  in  ray  parameters  (0.07  for  pP  and  0.12  for  sP).  The 
remaining  difficulty  is  in  estimating  (see  Figure  24).  If  we  assume  that  4t’,  we  obtain  the 
expected  value  of  -  3.6  Under  this  assumption,  which  is  supported  by  many  earthquake  studies, 
it  becomes  possible  to  determine  the  appropriate  combination  of  B  and  'k.  to  satisfy  the  pS 
waveforms  and  absolute  amplitude  levels.  The  average  amplitude  of  the  CANNIKIN  pS  phase 
from  21  observations  is  2083  mp  which  implies  that  ^=0.4  and  y.  -  7.3x  1 0  "cm^.  This  value  of 
y, agrees  remarkably  well  with  that  determined  by  Rayleigh  wave  modeling,  as  reported  by  Toksdz 
and  Kehrer  (1972)  who  obtained  y  .  -  7  2  v  1 0 "cm^  Thus,  we  set  the  absolute  level  of  CANNIKIN 
to  be  7.3  X  lO'^cm^  and  adjust  the  other  two  events  accordingly  to  obtain  the  RDP’s  given  in 
Figure  2lc.  If  we  follow  the  basic  scaling  arguments  of  Mueller  and  Murphy  (1971),  namely  that 
yield  (F)  is  proportional  to  effective  cavity  volume,  we  expect 

y.  =  (1) 

and 

A' -  (2) 
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Figure  23.  Comparison  of  observed  long  period  pS  waveforms  and  synthetics  using  different  values 
of  for  event  CANNIKIN.  The  source  model  and  velocity  structure  are  obtained  from  near  field 
modeling. 
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Figure  24.  Synthetic  pS  amplitudes  for  event  CANNIKIN,  where  the  first  peak  amplitude  is  plotted 
against  /*  for  different  values  of  B  and  v. 


The  two  scaling  constants,  Cj  and  C2,  can  then  be  estimated  by  regression  assuming  the  yields  of 
LOnIGSHOT,  MILROW  and  CANNIKIN.  The  results  are 

£-1=12.3  X  10* 

C2=4.7 

The  parameter  B  now  depends  on  depth  and  is  given  by 

Log  B  =  -0.05  X  1.3  Log  h 

These  formulas  prove  quite  effective  in  explaining  the  GASBUGGY-RULISON  tests  as  discussed 
later. 

Scaling  Relationships  gl  NTS:  The  data  set  at  NTS  is  much  larger  than  that  at  Amchitka  but  again 
consists  mostly  of  a  collection  of  strong  motion  and  teleseismic  seismograms.  There  also  exists  an 
abundance  of  free-field  data  at  distances  of  less  than  a  few  kilometers.  Most  of  this  type  of  data 
is  for  small  events  whose  yield  is  less  than  15  kilotons.  The  few  measurements  that  do  exist  for 
larger  events  are  inconsistent  in  most  situations  (Murphy,  1991). 

The  strong  motion  records  at  distances  beyond  a  few  kilometers  are  much  more  complicated 
at  NTS  than  at  Amchitka,  presumably  because  of  the  complex  local  geology  (Barker  el  al.,  1991). 
However,  the  initial  motions  of  the  type  discussed  earlier  for  CANNIKIN  (Figure  20)  can  still  be 
modeled  using  the  structure  derived  in  Barker  el  al.  Figure  25a  through  25e  display  complete 
seismograms  for  the  best  recorded  events.  The  columns  on  the  left  are  appropriate  for  the  5=1 
assumption  used  earlier  in  the  extensive  intercorrelation  exercise  (Lay  ct  al.,  1984).  The  columns 
on  the  right  correspond  to  using  a  slightly  different  scaling  law  than  Barker  ei  al.  and  allowing 
for  an  adjustment  in  B  to  fit  the  estimated  long-period  source  strengths.  The  synthetics  generated 
in  Figure  25  were  produced  by  Filon-AS,  a  frequencx -wavenumber  code.  The  crustal  model  is 
given  in  Table  1. 
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Figure  25  (a).  Near-in  observations  (top  bold  traces)  and  predictions  (bottom  traces)  for  event 
BOXCAR.  The  station  names  and  distances  are  printed  on  the  left;  amplitudes  are  printed  on  the 
right,  and  values  of  K,  V.(PSI),  and  B  are  printed  above  each  column. 
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Figure  25  (b).  Near-in  observations  (top  bold  traces)  and  predictions  (bottom  traces)  for  event 
MAST.  The  station  names  and  distances  are  printed  on  the  left;  amplitudes  are  printed  on  the 
right,  and  values  of  K,  V.(PSI),  and  B  are  printed  above  each  column. 
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Figure  25  (d).  Near-in  observations  (top  bold  traces)  and  predictions  (bottom  traces)  for  event 
Halfbeak.  The  station  names  and  distances  are  printed  on  the  left;  amplitudes  are  printed  on 
the  right,  and  values  of  K,  H'.(PSI),  and  B  are  printed  above  each  column. 
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Figure  25  (e).  Near-in  observations  (top  bold  traces)  and  predictions  (bottom  traces)  for  event 
SCOTCH.  The  station  names  and  distances  are  printed  on  the  left;  amplitudes  are  printed  on  the 
right,  and  values  of  K,  y.(PSI),  and  B  are  printed  above  each  column. 
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TABLE  1.  Velocity  Structure  Model 


Layer 

a,  km/s 

(3,  km/s 

p,  g/cm* 

Thickness,  km 

1 

2.30 

1.35 

1.90 

0.36 

2 

2.80 

1.50 

2.00 

0.40 

3 

3.30 

1.50 

2.25 

0.70 

4 

4.00 

1.90 

2.30 

0.70 

5 

4.60 

2.00 

2.40 

0.75 

6 

5.30 

2.50 

2.50 

0.80 

7 

5.50 

2.95 

2.70 

2.25 

8 

6.10 

3.50 

3.00 

10.00 

9 

7.00 

4.00 

3.00 

10.00 

Model  from  Hartzell  et  al.  (1983) 


The  long-period  V.  estimates  are  less  certain  for  NTS  events  for  several  reasons.  First,  the 
F  factor  or  tectonic  release  is  much  higher  than  at  Amchitka,  making  it  more  difficult  to  correct 
the  surface  waves  as  well  as  the  long-period  SV- waves  for  contamination.  Secondly,  the  geologic 
structure  is  more  complex  in  the  source  region,  so  that  individual  shots  sample  different  source 
parameters  (Murphy,  1989).  Results  from  Rayleigh  wave  modeling  by  Stevens  (1986)  is  given  in 
Table  2  along  with  V.  estimates  from  the  intercorrelation  technique  (Lay  et  al.,  1984a).  Note  that 
V. estimated  from  wj,  is  consistently  smaller  than  '•'.estimated  from  M^.  The  ratio  is  about  three 
to  one  but  with  considerable  scatter.  The  depth  effect  on  B  is  not  particularly  obvious  in  this  data 
set.  Note  that  the  event  CHESHIRE  is  deeper  than  ESTUARY  although  its  LP/SP  ratio  is  smaller. 
This  does  not  support  the  B  vs.  depth  dependence  discussed  earlier  in  connection  with  Amchitka 
data.  The  Rayleigh  wave  strength  given  by  Stevens  (1986)  is  in  general  agreement  with  those  given 
by  Given  and  Mellman  (1986)  such  that  a  factor  of  about  three  appears  appropriate  in  the  '•'. 


TABLE  2.  Strength  Estimates  from  Surface  Waves  and  Body  Waves 


Events 

Date 

'I'-f^b) 

Depth  (km) 

Ratio 

Scotch 

05/23/67 

0.91 

0.19 

0.98 

4.8 

Stinger 

03/22/68 

0.48 

0.15 

0.67 

3.2 

Sled 

08/29/68 

0.53 

0.30 

0.73 

1.8 

Almendro 

06/06/73 

1.9 

0.68 

1.06 

2.8 

Tybo 

05/14/75 

1.0 

0.46 

0.77 

2.2 

Stilton 

06/03/75 

0.48 

0.24 

0.73 

2.0 

Mast 

06/19/75 

1.2 

0.49 

0.91 

2.4 

Cheshire 

02/14/76 

2.0 

0.35 

1.17 

5.7 

Estuary 

03/09/76 

2.0 

0.35 

0.87 

5.7 

Pool 

03/17/76 

0.43 

0.44 

0.88 

1.0 

Following  the  results  of  strong-motion  modeling  assuming  B=\  (Hartzell  et  al.,  1983)  we 
used  the  parameterizations  given  earlier  in  equations  (I)  and  (2)  obtaining 

Cl  «  20.1  X  10* 

C2  =  4 

where  Ci  has  been  increased  by  a  factor  of  3.  These  formulas  apply  when  using  B*l.  If  we  used 
the  corner  frequency  scaling  relative  to  Amchitka  we  obtain  the  synthetic  comparisons  given  in 
Figure  26.  These  begin  to  produce  too  much  short-period  energy  at  the  nearest  stations  and  do 
not  fit  the  acceleration  data  (Helmberger  el  al.,  19913).  Thus,  the  corner  frequency  difference 
between  the  two  sites  appears  real.  We  now  have  MILROW  and  BOXCAR  with  K’s  of  9  and  7 
and  H'.’s  of  1.8  X  lO'^cm*  and  3.6  x  I0“cm*.  The  difference  in  S',  accounts  for  a  Mfs=0.48  which 
we  essentially  constructed. 

Off-Test  Site  Events:  The  local  and  teleseismic  data  from  various  off-test  site  events  have  been 
extensively  studied  in  attempts  to  explain  the  scatter  from  standard  curves.  For  example, 

LRSM  stations  were  established  at  FAULTLESS,  RULISON,  and  GASBUGGY  epicenters  to  measure 
relative  attenuation  to  explain  6mb’s  with  respect  to  6f‘  differences.  These  results  generally  do  not 
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Figure  26  (a).  Near-in  observations  (top  bold  traces)  and  predictions  (bottom  traces)  for  event 
BOXCAR  using  Amchitka  scaling.  The  station  names  and  distances  are  printed  on  the  left; 
amplitudes  are  printed  on  the  right,  and  values  of  K,  V.  (PSI),  and  B  are  printed  above  each 
column. 
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Figure  26  (b).  Near- in  observations  (top  bold  traces)  and  predictions  (bottom  traces)  for  event 
MAST  using  Amchitka  scaling.  The  station  names  and  distances  are  printed  on  the  left;  amplitudes 
are  printed  on  the  right,  and  values  of  K,  V.  (PSI),  and  B  are  printed  above  each  column. 
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Figure  26  (c).  Near-in  observations  (top  bold  traces)  and  predictions  (bottom  traces)  for  event 
INLET  using  Amchitka  scaling.  The  station  names  and  distances  are  printed  on  the  left;  amplitudes 
are  printed  on  the  right,  and  values  of  K,  'k.(PSI),  and  B  are  printed  above  each  column. 
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Figure  26  (d).  Near-in  observations  (top  bold  traces)  and  predictions  (bottom  traces)  for  event 
HALFBEAK  using  Amchitka  scaling.  The  station  names  and  distances  are  printed  on  the  left; 
amplitudes  are  printed  on  the  right,  and  values  of  K,  V.  (PSI),  and  B  are  printed  above  each 
column. 
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Figure  26  (e).  Near-in  observations  (top  bold  traces)  and  predictions  (bottom  traces)  for  event 
SCOTCH  using  Amchitka  scaling.  The  station  names  and  distances  are  printed  on  the  left; 
amplitudes  are  printed  on  the  right,  and  values 
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appear  to  explain  the  anomalies  (Der  et  al.,  1980).  For  example,  the  attenuation  at  the  RULISON 
site  is  nearly  the  same  as  at  GASBUGGY,  and  attenuation  at  FAULTLESS  is  the  same  as  it  is  at 
NTS.  Thus,  anomalies  in  remain  unexplained.  On  the  other  hand,  Afg  values  prove  effective 
in  estimating  these  yields  as  discussed  by  Yacoub  (19*3). 

A)  RULISON  and  GASBUGGY.  a  Test  Qf  Scaling:  The  two  non-NTS  nuclear  blasts,  RULISO.N 
and  GASBUGGY,  provide  a  particularly  useful  test  case  for  determining  the  transportability  of 
different  methods  of  estimating  yields  (see  Table  3).  These  two  nuclear  events  were  detonated 
304  km  apart  in  nearby  sedimentary  basins  in  New  Mexico  and  Colorado.  Although  the  announced 
yield  of  RULISON  was  40  kt  and  that  of  GASBUGGY  was  29  kt,  GASBUGGY  has  a  larger  mt- 
The  most  thorough  analysis  of  this  discrepancy  was  presented  by  Murphy  and  Archambeau  (1986). 
Although  the  amplitude  differences  between  the  two  events  are  of  the  order  that  might  be  expected 
from  a  difference  in  /’  of  0.4  sec  between  the  two  events,  their  analysis  of  the  P-wave  displacement 
spectra  of  the  two  events  indicated  that  the  effective  /^for  the  two  events  are  essentially  the  same, 
in  agreement  with  Der’s  (1980)  results.  Murphy  and  Archambeau  present  evidence  that  the 
RULISON-GASBUGGY  anomaly  is  principally  caused  by  tectonic  release  associated  with 
RULISON.  They  suggest  that  this  release  of  tectonic  energy  was  oriented  in  such  a  fashion  as  to 
destructively  interfere  with,  and  reduce  the  amplitudes  of,  the  teleseismic  short-period  P-waves 
for  RULISON.  The  magnitude  of  their  proposed  release  is  such  that  the  amplitude  of  the  tectonic 
P-waves  is  roughly  one  half  of  the  amplitude  of  the  explosion  P-waves. 

TABLE  3.  Explosion  Source  Parameters 


Event 

Location 

Date 

Yield  (kt) 

Depth  (m) 

GASBUGGY 

36.68°  N,  107.21°  W. 

10/12/1967 

29 

1292 

RULISON 

39.41°  N.  107.95°  W. 

10/09/1969 

40 

2573 
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Although  the  tectonic  release  hypothesis  cannot  be  ruled  out,  the  amount  of  tectonic  release 
rec  uired  to  produce  this  much  interference  is  difficult  to  support  from  the  results  presented  in 
Burdick  el  al.  (1991a).  Another  interpretation  is  that  this  reduction  in  RULlSON’s  amplitude  is 
caused  by  a  decrease  in  the  B  factor  as  suggested  by  the  above  scaling  relations  for  Amchitka. 

Using  the  Amchitka  scaling  rela*ions,  the  predicted  values  for  GASBUGGY  (//=1.29  km, 
y=29  kt)  are  B-0.64,  K«=31,  and  Y.  -  4.5x  10*cm^  or  Mo=2.6  x  IO22  dyne-cm.  For  RULISON 
(h=2.5~>  km,  y=40  kt),  the  predicted  values  are  5=0.27,  K=37,  and  Y.  -  5.9xr  lO’cm^or  Mo=3.7  x 
10^2  dyne-cm  (Figure  27).  Synthetic  seismograms  calculated  using  these  source  functions  parameters 
and  the  crustal  structures  given  by  Murphy  and  Archambeau  predict  that  the  for  GASBUGGY 
is  0.23  units  higher  than  that  for  RULISON.  This  difference  is  primarily  caused  by  the  small  B 
value  for  RULISON. 

Although  the  Amchitka  scaling  laws  do  a  good  job  at  predicting  the  values  of  GASBUGGY 
and  RULISON  an  examination  of  the  recorded  RDP’s  of  GASBUGGY  (Figure  28)  shows  that  the 
prediction  is  not  perfect.  Only  one  of  the  RDP’s  shows  a  significant  overshoot.  The  mean  overshoot 
is  15%,  corresponding  to  a  5  of  0.32.  A  K  value  approximately  equal  to  27  with  a  V .  -  6.4.y  lO’cni^, 
were  needed  to  model  the  rise  time  (see  the  bottom  trace).  These  observations  of  the  direct  RDP 
are  rather  unique  in  that  gas-wells  were  available  for  subsurface  observations  (Perret,  1982).  They 
were  all  taken  at  a  slant  range  of  about  0.5  km,  well  outside  the  nonlinear  zone.  The  average  RDP 
at  the  bottom  of  Figure  28  shows  a  higher  y.and  a  lower  overshoot  than  predicted  by  any  scaling 
laws. 

The  moments  obtained  from  modeling  the  surface  waves  (Burdick  el  al.,  1991a)  are  2.0  x 
1022  and  2.1  x  IO22  dyne-cm,  which  are  somewhat  low  compared  to  both  the  observed  RDP  and 
predicted  RDP.  These  discrepancies  are  probably  within  modeling  errors,  but  the  evidence  for 
undershoot  in  this  experiment  is  strong. 
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Figure  27.  Reduced  displacement  potentials  for  PNE’s  GASBUGGY  (5=0.64,  K=3l,  V,=4.5xl0«) 
and  RULISON  (B=0.27,  K=37,  y,;=5.9xI09). 


Figure  28.  Upper  four  traces  display  the  observed  RDP’s  obtained  from  the  GASBUGGY  explosion, 
after  Ferret  (1982).  Bottom  trace  displays  a  scaled  RDF. 
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Bi  FAULTLESS  and  ALMENDRO.  Further  Test  fif  Scaling:  Since  the  F-waveforms  from  Pahute 
Mesa  events  have  been  studied  at  great  length  by  Burdick  and  his  associates  it  is  natural  to  compare 
data  from  the  FAULTLESS  event  with  a  Pahute  Mesa  event  (ALMENDRO)  having  comparable 
yield  (see  Figure  29).  Both  blasts  were  detonated  in  similar  environments.  ALMENDRO 
was  located  roughly  150  km  to  the  north.  The  pairs  of  seismograms  are  from  the  same  teleseismic 
WWSSN  stations.  It  appears  that  the  FAULTLESS  seismograms  are  anomalously  long  period  for 
such  a  small  event.  This  suggests  a  higher  t’^  but  several  detailed  attenuation  experiments  indicate 
comparable  values  of  /^at  these  two  sites.  NTS  F- waveforms  for  the  large  events  such  as  BOXCAR 
are  similar  to  FAULTLESS  and  intercorrelation  of  FAULTLESS  using  BOXCAR  as  a  master  event 
produce  the  same  yields  (Lay  et  al.,  1983).  In  short,  the  F-waveforms  from  FAULTLESS  fall 
nicely  into  the  Pahute  Mesa  population,  but  are  too  strong.  Allowing  the  B  to  increase  by  roughly 
a  factor  of  two  for  FAULTLESS  can  explain  this  6mb  anomaly  as  discussed  in  Burdick  et  al.  (1991a). 
Since  the  A/g's  for  these  events  are  roughly  the  same,  we  consider  this  a  reasonable  explanation. 
Discussion:  In  previous  sections  we  have  reviewed  the  data  and  scaling  law  applications  for  various 
U.S.  shots.  In  this  section  we  discuss  the  yield  equations  put  forward  by  other  researchers  for  sites 
outside  US  and  how  we  think  these  should  be  altered  based  on  this  study.  Because  of  strong 
evidence  for  overshoot  variations  from  site  to  site  we  believe  that  magnitudes  based  upon  long-period 
excitation  are  the  most  reliable.  Thus,  we  begin  with  yield  estimation  based  on  A/g-  From  Sykes 
and  Wiggins  (1986)  we  obtain  the  Novaya  Zemlya  scaling  relation 

5 -0.971  log  ,o>  *2.16 

This  equation,  in  fact,  matches  the  Amchitka  data  better  than  does  the  Amchitka  scaling  relation 
presented  earlier, 

Al 5  "  log  loF  *  2.0 
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Figure  29.  A  comparison  of  seismograms  of  observed  WWSSN  F-waveforms  for  the  two  explosions 
ALMENDRO  (A)  and  FAULTLESS  (F)  as  a  function  of  azimuth  (same  station).  The  numbers 
indicate  the  ratios  of  peak  amplitudes  (FAULTLESS/ALMENDRO). 


when  the  most  accurate  yield  values  are  inserted.  This  does  not  greatly  change  the  modeling 
exercise  discussed  earlier  except  to  make  the  dependence  of  B  upon  depth  even  stronger.  The 
formula  by  Sykes  and  Cifuentes  (1984)  for  Eastern  Kazakh  is  nearly  the  same  as  above,  namely, 

Af5-0.951og  ,oy  +  2.16 

The  results  for  NTS  show  an  obvious  off-set  from  Amchitka  and  other  sites  as  discussed  earlier. 
If  we  compare  the  10  largest  events  fired  at  Shagan  River  versus  the  10  largest  events  fired  at  NTS 
after  the  1976  150  kt  threshold  (as  suggested  by  Figure  16  of  Stevens,  1986),  we  obtain 

Afl(NTS)-'XiU(ShaganRiver)-  0.63 

with  an  average  Mq  of  7.58  x  10^2  dyne-cm  at  Shagan  River  vs.  an  average  Mo  of  1.74  x  10” 
dyne-cm  at  NTS.  Our  formula  predicts  a  for  NTS  which  is  0.39  units  higher  than  that  for 
Amchitka.  This  suggests  that  the  H'.  scaling  laws  for  Amchitka  are  applicable  to  Soviet  tests  (as 
assumed  by  Sykes  and  his  colleagues).  Thus,  NTS  appears  to  have  a  unique  scaling  law  characterized 
by  undershoot  compared  to  any  other  site.  If  we  suppose  that  most  of  the  Novaya  Zemlya  tunnel 
events  are  fired  at  similar  depths,  we  conclude  that  the  smallest  events  are  over-buried  while  the 
larger  events  are  under-buried.  Thus,  at  Novaya  Zemlya  we  would  expect  small  B's  for  small 
events  and  large  B's  for  large  events,  or  mj,  increasing  with  increasing  yield  more  rapidly  than  at 
Amchitka  (see  Figure  30).  For  example,  suppose  we  compare  the  larger  events  at  Amchitka  (where 
the  data  is  the  most  complete)  with  comparable  events  at  Novaya  Zemlya.  For  the  largest  Novaya 
Zemlya  event  (October  27,  1973),  the  Sykes  formula  predicts  a  yield  of  2840  kt  for  the  observed 
\fs  of  5.5.  The  'k.  obtained  for  this  event  by  Cohee  and  Lay  (1988)  is  3.8  x  10“  cm’  (assuming 
the  SL8  attenuation  model  proposed  by  Anderson  and  Hart,  1978).  This  estimate  was  made  using 
a  /^of  4.8,  which  is  probably  high.  The  trade-off  between  N^.and  fji  was  shown  earlier  in  Figure 
24.  Reducing  tl  to  0.5  and  to  2.0  suggests  a  H'.  near  2  x  10“  cm’,  a  value  midway  between  that 
of  MILROW  and  CANNIKIN.  The  mb  for  the  Novaya  Zemlya  event  is  6.9  or  that  of  CANNIKIN. 
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Variations  in  t* 


Figure  30.  (Top  panel)  Effect  of  varying  on  events  LONGSHOT,  MILROW.  and  CANNIKIN 
while  holding  source  function  and  moment  constant.  The  solid  line  represents  actual  vs. 
values.  (Bottom  panel)  Effect  of  varying  B  for  same  events  while  holding  f*  and  moment  constant. 
The  crosses  represent  actual  vs.  Afg  values. 
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However,  allowing  B  to  grow  for  this  large  Novaya  Zemlya  event  easily  reduces  its  y.  to  a  value 
near  that  of  MILROW.  In  short,  it  is  very  difficult  to  use  the  mb’s  to  directly  estimate  yields 
while  dealing  with  these  trade-offs  of  /*  and  B. 

The  mb  measurement  thus  is  shown  to  have  several  problems  in  estimating  the  yield  of  bombs 
in  geological  regimes  which  do  not  have  numerous  test  shots  available.  These  include  unknown 
factors  of  corner  frequency,  overshoot,  tectonic  release,  and  /*.  Using  a  regional  Mg  measurement, 
when  possible,  should  serve  to  minimize  the  problems  caused  by  these  factors.  Such  factors  as 
explosion  coupling  will  still  remain  a  problem. 

In  most  instances,  however,  Soviet  mb’s  will  still  have  to  be  used.  Perhaps  the  simplest  use 
of  the  Soviet  mb’s  is  to  directly  relate  them  to  Soviet  A/g  vs.  yield  relations.  By  taking 

m,-  1.05Af5+  1.63 

and  substituting 

Af3-0.971og,o>'*2.16 

we  obtain 

1 .021  log,o>'  +  3.9. 

Thus  for  events  that  have  high  T-factors  we  can  still  estimate  the  yields  from  mb.  This  equation 
estimates  the  mb  of  a  150  kt  event  at  Shagan  River  to  be  6.1  as  observed  during  the  recent  PVE. 
However,  instead  of  adjusting  the  U.S.  curve  by  applying  the  6mb  bias  correction,  we  adopt  the 
scaling  relation  from  that  of  Amchitka  and  explain  the  difference  in  mb  by  overshoot  and  attenuation. 
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Q  STUDIFS  PART  1  -  EFFECT  OF  Q  ON  THE  REGIONAL  Pn  and_SiL  COMPOSITION 

RECORDED  IN  U.S.  ANP  U^.S.R 

Introduction:  With  the  installation  of  the  broad-band,  high  dynamic  range  IRIS  instruments,  it 
has  become  possible  to  compare  the  regional  waveforms  of  earthquakes  and  explosions  at  magnitudes 
3  to  6.  The  dynamic  range  feature  of  these  systems  allows  the  comparison  of  the  relatively  waek 
body  waves  (Pnl)  with  the  stronger  surface  waves.  The  broad-band  feature  allows  the  examination 
of  the  frequency  content  of  particular  phases  (Pn  and  Sn)  to  address  the  Q-issue.  Computational 
methodologies  have  advanced  in  recent  time  which  together  with  the  development  of  recent  fast 
computational  facilities  has  made  it  feasible  to  investigate  the  response  of  a  laterally  varying  crustal 
medium  within  a  reasonable  time  frame.  Among  other  methods,  the  generalized  ray  theory  is  a 
most  widely  used  method  for  analysis  of  the  composition  of  seismic  phases  that  the  source  process 
and  the  propagation  effect  together  make  up  at  a  receiving  station.  While  this  method  is  applied 
to  a  flat-layered  crustal  structure  in  most  studies,  the  method  has  recently  been  extended  to  include 
range-dependent  structure  (Helmberger,  personal  communication).  It  is  a  computationally  fast 
technique,  but  requires  many  generalized  rays  to  be  tracked  between  the  source  and  the  receiver 
for  a  regional  waveguide.  As  the  crustal  medium  becomes  complicated,  the  method  can  quickly 
become  quite  cumbersome  with  the  process  of  just  tracking  the  rays.  The  anelasticity  of  the 
waveguide  is  applied  to  resulting  response  of  many  generalized  rays  in  terms  of  t’  which  is  an 
average  estimate  of  seismic  attenuation.  In  actuality,  seismic  waves  attenuate  in  different  amounts 
depending  on  the  material  property  and  should  be  treated  as  an  intrinsic  behavior  of  the  medium. 
Thus,  Q  (quality  factor)  should  be  defined  for  each  crustal  layer  for  both  compressional  and  shear 
waves. 

At  high  frequencies  the  effects  of  scattering  in  the  crust  become  so  intense  that  only  statistical 
properties  of  waveforms  are  meaningful.  The  receiver  related  crustal  structure  can  be  complex 
and  may  cause  additional  complexity  beyond  that  already  caused  by  the  triplications  of  seismic 
waves  due  to  the  gradual  velocity  increase  of  the  crust-mantle  transition  zone.  The  presence  of  a 


gradient  structure  near  the  receiver  changes  the  timing  of  surface  reflected  phases  and  their 
reflection  coefficients.  Consequently,  the  waveforms  become  complex  and  generalized  ray  theory 
can  be  used  to  identify  the  significant  arrivals  within  the  composition  of  short-period  signals  (Saikia 
and  Burdick,  1991 The  most  suitable  computational  method  to  apply  for  investigating  the  effect 
ul'  u  crustal  waveguide  of  this  nature  is  the  method  of  frequency-wavenumber 
integration/reflectivity.  It  allows  evaluation  of  full  medium  response  where  the  intrinsic  aspects 
of  the  attenuation  of  the  medium  can  be  specified  and  effects  on  the  Pn  and  Sn  composition  can 
be  investigated. 

In  a  recent  study,  Saikia  and  Burdick  (1991)  showed  that  short-period  Pni  waves  (period  as 
short  as  2  s)  are  stable  and  can  be  modeled.  They  studied  many  observations  from  Nevada  Test 
Site  (NTS)  explosions  recorded  at  regional  distances  of  200  to  420  km  and  modeled  the  P„i 
waveforms  using  a  deterministic  crustal  waveguide.  The  sources  of  these  waveforms  were  shallow. 
Also,  the  sources  were  predominantly  isotropic,  and  the  portion  of  P„i  waves  which  was  included 
in  the  Pn  and  Pg  waves  had  a  duration  of  about  30  s  and  was  dominated  by  compressional  waves. 
To  understand  the  observed  data,  they  used  the  frequency-wavenumber  algorithm  to  compote  the 
explosion  generated  Pni  waves  for  several  canonical  crustal  models  and  selected  a  crustal  model 
based  on  the  agreement  between  the  data  and  the  synthetic  seismograms.  The  method  was  then 
utilized  to  understand  the  composition  of  the  Pg  wave  group  which  was  constituted  of  phases  like 
PmP,  pPmP,  2PmP,  PmS,  pPmS,  PmPSmP,  PmSPmP  etc  for  realistic  models.  Pg  is  a  wave  group 
whose  frequency  content  is  widely  used  to  discriminate  events.  In  this  study,  we  have  taken  a 
similar  strategy  to  investigate  the  broadband  composition  of  Pnj  and  Sni  seismograms  recorded  in 
the  North  American  continent  at  regional  distance  from  double-couple  sources  and  of  the  P  and 
S  waveforms  that  are  recorded  within  the  Soviet  Union.  We  shall  mainly  focus  on  identifying  the 
rays  important  to  model  the  regional  waves  within  the  S  wave  window  and  investigate  how'  the 
intrinsic  Q  model  affects  the  waveform  composition. 


62 


Data:  For  the  U.S.  study,  we  used  a  set  of  three-component  broadband  seismograms  recorded  at 
Har/ard  (HRV)  station  at  a  distance  of  640  km  from  the  Saguenay  earthquake  of  November  25, 
1988  (Figure  31a).  These  seismograms  were  recorded  on  a  Streckeisen  seismometer.  We  selected 
these  seismograms  because  many  features  recorded  on  the  seismograms  were  successfully  modeled 
by  Zhao  and  Helmberger  (1991).  Beginning  with  their  crustal  model,  we  have  directed  our  study 
towards  the  modeling  the  high-frequency  details  observed  in  the  waves  and  the  composition 
of  waves  identified  as  and  sSn  by  Zhao  and  Helmberger  (1991)  using  a  multiple  source  model. 
A  similar  study  was  directed  towards  the  modeling  of  the  broadband  seismograms  recorded  within 
the  Soviet  Union.  We  selected  a  set  of  three -component  seismograms  recorded  at  GARM  from  an 
earthquake  which  originated  at  a  distance  of  200  km  at  an  azimuth  of  290°  on  May  4,  1989  (38.73°N 
and  78.5°W,  Figure  31b).  Unlike  for  North  American  earthquakes,  the  waveforms  of  very  few 
Soviet  Union  earthquakes  have  been  modeled.  Thus,  it  is  necessary  to  develop  a  starting  crustal 
model  even  to  obtain  a  first-order  agreement  between  data  and  synthetic. 

Modeling  ftf  HRV  Seismograms  from  the  Saguenay  Earthquake:  Figure  32  shows  the  broadband 
displacements  recorded  at  Harvard  station.  To  investigate  the  influence  of  crustal  structure  on  the 
various  significant  phases  of  the  window,  we  have  started  with  the  crustal  model  shown  in 
Figure  33  by  dotted  lines.  This  model  extends  from  the  surface  to  a  half  space  at  a  depth  of  55 
km.  The  major  velocity  discontinuity  is  at  a  depth  of  35  km  where  the  P  velocity  jumps  from 
6.71  km/s  to  8.1  km/s  and  the  S  velocity  jumps  from  3.82  km/s  to  4.7  km/s.  Zhao  and  Helmberger 
(1991)  used  a  reflectivity  code  (Mallick  and  Frazer,  1988)  to  compute  the  n.edium  response  and 
used  an  elastic  crustal  structure  to  model  the  data.  They  used  a  Qp (shear- wave  quality  factor)  of 
6200  and  stated  that  a  lower  value  of  Qgis  not  required  to  match  the  recorded  wave  form,  although 
the  conventional  wisdom  is  that  for  eastern  North  America  Q@is  of  the  order  of  300  (Hwang  and 
Mitchell,  1987).  Based  on  this  published  information,  we  started  to  look  for  certain  phases  within 
the  Pni  regime  for  which  the  agreement  between  the  data  and  synthetic  can  be  improved  and  in 
the  process  to  learn  more  about  the  regional  waveguide. 
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Figure  31.  Geographical  location  of  (a)  HRV  (HarvardO  station  (solid  triangle)  and  November  25, 
1988  Saguenay  earthquake  (solid  star)  and  (b)  GAR  (Garm)  station  (solid  triangle)  and  May  4,  1989 
USSR  earthquake  (solid  star). 


Brccdbond  Displacement  recorded  of  Horvard 
Stoficn  from  1988—11-25.  Soguenoy  Earthquake 
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Figure  32.  Broadband  three-component  displacement  seismograms  as  recorded  by  Harvard  station 
from  the  1988,  November  25  Saguenay  earthquake.  The  original  seismograms  were  integrated. 
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ENA  Regional  Velocity  Models 
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Figure  33.  Regional  crustal  model  developed  by  modeling  the  broad-band  seismograms  recorded 
at  Harvard  station  from  the  November  25,  1988  Saguenay  earthquake.  The  final  model  is  shown 
by  the  solid  lines.  The  model  shown  by  the  dotted  lines  is  the  initial  crustal  model  developed  by 
Zhao  and  Helmberger,  1990. 
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The  phases  marked  as  S„  and  sS„  show  the  greatest  misfit  between  the  data  and  the  synthetics 
conputed  by  Zhao  and  Helmberger  (1991)  (see  their  Figure  16).  The  synthetic  seismograms  are 
definitely  of  lower  frequency.  So  our  initial  attempt  was  to  understand  what  part  of  the  crustal 
waveguide  would  be  most  critical  in  development  of  these  waveforms.  In  the  present  calculation, 
we  used  the  frequency-wavenumber  integration  method  and  set  the  nyquist  frequency  at  10  hz. 
We  computed  theoretical  seismograms  for  eight  fundamental  faults  and  used  a  focal  mechanism 
with  a  dip  65®,  a  rake  of  78°  and  a  strike  of  323°  to  predict  the  vertical,  radial  and  tangential 
component  seismograms.  These  synthetics  were  used  to  compute  both  the  point  and  multiple  source 
seismograms  and  the  corresponding  vertical  component  seismograms  are  shown  in  Figure  34.  The 
source  model  contained  three  sub-sources,  with  seismic  moments  of  1.55x1024,  1.45x1024  and 
1.95x1 024  dyne-cm  respectively!  (after  Zhao  and  Helmberger,  1991).  The  second  source  was  delayed 
by  0.65  s  and  the  third  source  by  1.45  s  from  the  first  source  to  account  for  the  propagation  of 
the  rupture  front.  The  first  source  was  represented  with  a  source  time  function  defined  by  a 
trapezoid  of  0.4s  rise  time,  0.05s  of  follow-on  time  and  Q,25s  of  healing  time.  Similarly,  the  second 
and  third  sources  were  convolved  with  trapezoids  of  (0.2s,  0.15s,  0.15s)  and  (0.1s,  0.3s,  0.2s), 
respectively.  We  also  show  the  synthetic  seismograms  generated  by  Zhao  and  Helmberger  (1991) 
in  Figure  35  using  a  nyquist  frequency  of  4  hz  so  that  a  direct  comparison  can  be  made  with  those 
shown  in  Figure  34.  The  frequency  content  in  the  Pnl  waves  of  these  seismograms  is  not  as  rich 
as  those  Pnl  waves  shown  in  Figure  34  where  the  high  frequencies  are  the  result  of  derived  source 
complexity. 

In  Figure  36,  we  compare  the  vertical  and  radial  component  showing  just  the  Pnl  portion  of 
the  seismograms  computed  using  the  parameters  of  multiple  sources.  The  high-frequency  signals 
are  adequately  predicted  with  respect  to  those  observed  on  the  recorded  data.  The  seismograms 
computed  using  the  response  up  to  4  Hz  were  essentially  identical  to  these  seismograms. 
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Point— Source  Vs  Multiple— Source  Synthesis  of  the 
Harvard  Record  for  Vorious  Crustal  Models 
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Figure  34.  Comparison  Between  two  sets  of  displacement  seismograms  synthesized  using  point  and 
multiple  sources,  (a)  Point-source  displacements  -  the  upper  seismogram  is  computed  using  the 
model  response  of  Zhao  and  Helmberger  (1991)  and  the  bottom  seismogram  is  computed  using  the 
model  response  of  the  present  crustal  model,  and  (b)  multiple-source  displacements  for  the  two 
crustal  models. 
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Figure  35.  Comparison  between  data  and  synthetic  displacements  with  a  nyquist  of  4  Hz.  The 
seismograms  for  625  km  was  computed  using  a  different  velocity  crustal  model  (Figure  taken  from 
Zhao  and  Helmberger,  1991). 
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Figure  36.  Comparison  between  data  and  synthetic  displacements  with  a  model  response  up  to  a 
nyquist  of  10  Hz.  Note  the  developement  of  high  frequencies  and  agreement  between  phases 
marked  by  the  arrows. 
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Rav  Analysis  fi£  E|il  Seismograms:  in  this  section,  we  discuss  our  investigation  of  the  constituent 
phases  of  the  recorded  seismogram  at  Harvard  station.  The  basic  idea  is  to  investigate  the 
interaction  of  individual  ray  groups  in  creating  the  total  seismogram.  We  computed  generalized 
ray  seismograms  using  the  source  process  of  the  Saguenay  earthquake  for  several  groups  of  generalized 
rays.  In  Figure  37,  we  display  \ ertical-component  seismograms  of  these  ray  groups.  The  top  si.\ 
seismograms  are  normalized  to  their  maximum  amplitude.  All  the  PmP  and  SmS  rays  were  allowed 
to  reflect  from  each  interface  beneath  the  crust-mantle  boundary  including  the  reflection  from 
the  Moho  discontinuity.  The  total  response  of  these  PmP  and  SmS  rays  is  plotted  in  the  first 
seismogram.  The  geometric  arrivals  are  indicated  by  PmP  and  SmS  respectively.  The  S„  arrival 
is  small  and  is  preceded  by  a  refracted  phase  SP.  This  refracted  phase  had  developed  due  to  a 
critical  incidence  of  an  S  wave  on  an  interface  permitting  the  converted  P  phase  to  travel  along 
the  interface.  The  seismogram  in  the  second  row  is  for  the  sPmP,  a  ray  which  has  departed  from 
the  source  as  a  S  wave  and  then  converted  to  P  mode  at  the  free  surface.  The  amplitude  of  this 
ray  is  small.  The  next  seismogram  is  for  sSmS.  Both  the  geometric  and  head  waves  are  strong  for 
this  ray  group  and  contribute  significantly  to  the  total  seismogram.  The  next  two  seismograms  are 
for  the  SmSSmS  and  sSmSSmS  ray  groups.  Both  the  ray  groups  have  significant  contributions.  The 
sixth  seismogram  is  for  a  ray  group  identified  as  SmS'SmS.  The  rays  included  in  this  group  leave 
the  source  downward  and  reflect  from  each  interface.  The  reflections  are  turned  back  into  the 
lower  crust  again  at  the  Moho  discontinuity  before  they  are  reflected  back  to  the  receiver.  The 
contributions  from  these  rays  do  offer  a  significant  contribution  to  the  evolution  of  the  wave 
group.  The  seismogram  "Total"  is  the  result  of  direct  sum  of  the  upper  six  seismograms.  Having 
obtained  a  good  agreement  between  the  data  and  the  synthetics,  we  plotted  the  multiple-source 
frequency  wa'enumber  seismogram  computed  using  the  frequency- wavenumber  method  beneath 
the  total  response  for  a  direct  comparison.  This  comparison  produced  good  agreement  among  the 
dominant  features  within  the  so  called  waves". 
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Oistonce  -  R  =  640.0  Km  (Soguenoy  Epicenter  to  HARVARD) 


Figure  37.  Understanding  the  wavform  recorded  at  Harvard  station  using  the  ray  decomposition 
technique.  The  top  six  seismograms  are  for  the  individual  ray  groups.  The  seismogram  labelled 
"Total"  is  the  total  response  of  all  the  responses  of  upper  six  seismograms  and  the  comparison  with 
the  F-K  seismograms  shown  below  suggests  a  good  agreement  between  the  two  seismograms. 
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Thus,  we  have  extended  our  previous  study  (Saikia  and  Burdick,  1991)  on  the  deciphering 
of  he  ray  composition  of  Pg  waves  from  explosion  sources  to  earthquake  sources.  As  in  the  above 
study,  we  found  that  the  waveforms  within  the  Sn  group  can  be  studied  in  time  domain  in  terms 
of  a  basic  few  rays,  namely  the  SmS,  sSmS,  SmSSmS,  sSmSSmS  and  SmS'SmS  rays.  Since  these 
pha>es  leave  the  source  as  S  waves,  they  are  not  excited  by  the  explosion  source.  Therefore,  the 
only  phases  that  may  arrive  within  the  Sn  widow  from  a  pure  isotropic  source  are  the  P  waves  that 
are  converted  to  S  waves. 

Effect  of  Intrinsic  Q  gE  En  and  Sn:  It  is  expected  that  anelasticity  will  influence  the  frequency 
content  of  the  Pn  and  Sn  waves.  In  this  study,  we  wish  to  investigate  the  effects  on  regional  Pn 
and  Sn  composition  related  to  intrinsic  Q,  especially  the  changes  relative  to  the  seismograms  of  an 
elastic  medium.  Our  objective  is  to  determine  the  group  of  waves  that  is  most  sensitive  to  the 
variation  in  Q.  We  assume  that  Q  varies  significantly  within  the  upper  crustal  medium.  An  initial 
anelastic  calculation  was  performed  using  a  Qeof  250  in  the  upper  35  km  of  the  crust.  Qa  was 
assumed  to  be  twice  of  <?&  Below  35  km,  the  entire  medium  was  treated  as  elastic.  In  an  another 
Q  model,  Qg'vas  lowered  subtantially  from  250  to  50  within  the  upper  5  km  of  the  crust.  Figure 
38  shows  a  comparison  of  the  recorded  broadband  vertical  seismograms  at  Harvard  station  (top 
seismogram)  with  three  synthetic  seismograms  resulting  from  the  elastic  and  two  anelastic 
propagation  media.  The  second  seismogram  is  mutiple-source  seismogram  computed  using  the 
elastic  model.  While  the  frequency  content  of  the  signals  marked  by  arrows  compare  well  with 
the  data,  the  surface  waves  are  of  much  higher  frequency.  The  third  seismogram  is  synthesized 
using  the  anelastic  model  of  Qa  equal  to  250  throughout  the  35  km  thick  crust.  The  fourth 
seismogram  is  synthesized  with  a  Qnof  50  as  discussed.  The  bottom  two  seismograms  are  remarkably 
similar.  The  frequency  content  of  the  surface  waves  or  other  multiple  bounces  following  the  third 
arrow  are  in  better  agreement  with  the  data  although  the  fundamental-mode  Rayleigh  wave  appears 
longer  in  period. 
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In  Figure  38,  we  show  only  the  first  7Ss  of  the  seismograms  plotted  in  the  previous  figure. 
Clearly  the  similarities  between  data  and  synthetics  are  striking,  suggesting  that  the  waveshape  of 
the  P„|  waves  is  least  affected  by  the  variation  in  Q  model.  Figure  39  shows  a  comparison  of  both 
Pn  and  Sn  wave  groups  separately  as  a  function  of  different  Q  models.  The  seismograms  shown 
in  the  two  boxes  are  the  windowed  Pn  and  Sn  waves  shown  in  the  upper  two  seismograms.  Included 
in  the  Pn  panel  are  the  classical/turning  ray  Pn  and  PmP  wave,  and  in  the  Sn  panel  are  the 
classical/turning  ray  Sn,  surface  refracted  SP,  sSn  and  SmSSn  waves.  As  seen  in  the  previous 
figure,  the  P  waveforms  are  strikingly  similar  but  the  waveforms  drawn  in  the  Sn  box,  especially 
the  sSn  and  SmSSn  phases,  show  a  large  dependence  on  the  Q  model. 

This  study  pertains  to  the  Saguenay  earthquake  which  is  a  deep  event,  depth  27  km.  We 
would  expect  explosions  to  be  shallow.  But  the  Pn  generated  from  such  shallow  sources  spends 
only  a  slightly  longer  path  in  the  low  Q  materials.  Thus  it  can  be  assumed  that  the  initial  Pn  waves 
from  the  two  sources  will  be  affected  in  an  identical  manner  by  anelasticity.  However,  the  phases 
like  sSmS  and  SmSSmS  are  more  dramatically  affected  by  the  Q  model.  Therefore,  any  spectral 
ratios  of  Pn  to  SP+sSn;  pPn  to  SP+Sn  etc.,  may  provide  a  reasonable  measure  of  Q. 

Modeling  fif  GARM  Seismograms  from  Mav  1989  USSR  Earthquake;  In  this  study  we  have 
used  a  set  of  three-component  seismograms  recorded  at  the  Garm  station  from  an  earthquake  of 
May  4,  1989  (latitude;  39.436°N  and  longitude;  75.35°E,  h=35  km,  ISC).  The  station  is  located  at 
a  distance  of  200  km  from  the  source.  Figure  40  shows  the  recorded  displacements  processed  from 
the  broadband  velocity  seismograms.  A  high-pass  filter  was  applied  to  remove  the  long-period 
effects.  The  crustal  structure  encountered  by  the  wavefield  along  its  propagation  path  is  complex 
which  is  reflected  in  the  waveform.  To  begin  to  understand  the  waveforms,  it  was  necessary  to 
develop  a  crustal  structure. 

Our  strategy  for  developing  the  crustal  model  was  to  begin  with  the  tangential  component 
seismogram  because  of  the  simplicity  of  the  observed  displacement.  This  component  contains  only 
three  distinct  individual  arrivals  as  marked  by  the  arrows.  We  used  generalized  ray  theory  to 
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Effect  of  Crustol  Q  on  the  Whole  Regionol  Waveform 
ot  HARVARD  STATION  -  R  =  640  KM 
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Figure  38^  Comparison  between  broad-band  vertical  component  data  at  Harvard  station  and 
corresponding  synthetic  seismograms  computed  for  different  Q-models  including  the  elastic  model. 
Also  shown  is  the  first  75  seconds  of  each  seismogram  for  each  Q  model. 
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(o)  Da(o,  (b)  Elastic,  (c)  Crustol  0  =  250.  {d)  Upper  Crustol  (5km)  0=50 


Figure  39.  Comparison  between  data  broad-band  vertical  component  Pn  and  Sn  data  and 
corresponding  synthetic  seismograms  computed  for  different  Q-models  including  the  elastic  model. 
The  effect  of  Q  on  the  Sn  waves  is  more  pronounced  compared  to  the  effect  observed  on  the  Sn 
waves. 
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Broodbond  Displacement  recorded  ct  GARM  Stotion 
from  1989-05-04  Ecrthquoke  In  Soviet  Union;  AZ  =  292 
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Figure  40.  L'loadband  three-component  displacement  seismograms  as  recorded  by  Garm  station 
from  the  USSR  earthquake  of  May  4,  1989.  The  original  seismograms  were  integrated. 
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synthesize  this  component  for  various  crustal  models  and  several  source  depths.  The  best  prediction 
was  obtained  for  a  source  depth  of  25  km  for  the  preliminary  crustal  model  shown  in  Figure  41. 
We  used  the  following  focal  mechanism:  dip=75°,  slip=-135“  and  a  strike=32°.  The  structure 
contains  two  major  discontinuities  representing  the  Conrad  and  the  Moho.  We  succeeded  in  modeling 
the  tangential  displacement  using  only  three  arrivals,  the  direct  SH  and  two  reflections  from  the 
discontinuities.  The  frequency-wavenumber  seismograms  were  computed  with  these  parameters 
and  were  compared  with  the  recorded  data.  Figure  42  shows  the  comparison  between  the  data  and 
the  synthetic  seismograms  for  the  vertical  and  tangential  components.  We  are  successful  in  producing 
a  good  agreement  between  the  data  and  the  synthetic  for  the  tangential  motion.  We  mark  the 
individual  arrivals  in  the  synthetic  and  show  their  correspondence  with  the  data  by  the  thin  arrows. 
The  vertical  component  show  agreement  in  the  arrivals  of  Pn,  PmP,  ScS  and  SmS  phases.  The 
signal  bracketed  within  the  window  of  the  vertical-component  synthetic  seismogram  has  a  similar 
character  in  the  frequency  content  to  the  signal  bracketed  within  the  data  window.  Figure  43 
shows  a  generalized  ray  seismogram  using  the  three  rays.  The  vertical  synthetic  seismogram  has 
a  strong  SmS  which  is  smeared  out  in  the  data  due  to  the  interaction  with  the  physically  more 
complicated  crust  in  the  region.  The  vertical  component  of  the  recorded  seismogram  is  also 
dominated  by  long-period  P^i  signals  shown  by  the  solid  window.  These  waves  are  also  observed 
in  the  synthetic  seismogram  but  arriving  at  Garm  with  a  fast  velocity.  We  also  investigated  the 
effect  of  a  possib'e  linear  velocity  gradient  near  the  free  surface  to  determine  if  such  a  velocity 
distribution  would  account  for  the  mismatch  between  the  data  and  the  synthetics  within  bracketed 
the  window  shown  in  Figure  42.  We  discretized  the  top  ten  kilometers  of  the  crust  into  ten  layers 
of  equal  thickness  and  allowed  a  P-wave  velocity  increase  from  4.5  km/sec  to  5.5  km/sec  from 
the  surface.  The  S-wave  velocity  within  each  layer  had  a  ratio  of  1.73  to  the  P-wave  velocity. 
This  seemed  to  be  a  particularly  reasonable  explanation  for  the  small  complex  phases  between  the 
major  arrivals.  However,  the  synthetic  seismograms  computed  using  this  surface  gradient  did  not 
improve  the  fit  to  a  significant  d^'gree. 
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Figure  41.  Regional  velocity  used  for  modeling  GARM  seismogram. 


Preliminary  Modeling  of  Broadbond  Displacement 


Figure  42.  Comparison  is  shown  betwe'*n  the  data  and  the  synthetic  seismogram  for  the  vertical 
and  tangential  components.  The  signal  marked  by  arrow  l  is  a  phase  reflected  from  the  Conrad 
and  by  arrow  2  is  form  the  Moho  discontinuity.  The  phase  marked  by  arrow  3  on  the  tangential 
component  is  the  direct  arrival.  The  signal  within  the  window  shows  a  possible  correlation  between 
the  data  and  synthetics. 
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Figure  43.  A  vertical  component  seismogram  computed  using  direct  P  and  S,  PcP  ScS,  and  PmP 
and  SmS  phases.  These  arrivals  car.  distinctly  be  observed  on  the  recorded  seismograms. 
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Effect  oi  Intrinsic  Q  fifi.  Si  The  objective  of  this  study  is  to  determine  a  level  of  agreement  in 
the  frequency  content  between  the  recorded  and  the  predicted  SH  waves.  Similar  to  the  study  of 
Harvard  seismogram,  we  tried  two  Q  models  and  the  agreement  between  the  data  and  the  synthetic 
is  shown  in  Figure  44.  The  elastic  model  predicts  a  higher  frequency  content  for  both  the  direct 
and  conrad  reflected  SH  waves  (marked  by  the  arrows).  It  appears  <?b  lying  in  between  100  and 
250  would  predict  the  frequency  content  of  these  two  phases  in  better  agreement  with  the  data. 
The  Qpfor  the  Harvard  seismogram  was  close  to  250.  Thus,  Q  is  smaller  along  this  path  in  Garm. 
In  Figure  45,  we  show  four  vertical  component  seismograms  recorded  at  the  Garm  station  from 
four  separate  earthquakes  of  magnitude  ranging  from  4.7  to  5.0.  Two  of  the  seismograms  are 
recorded  at  an  azimuth  of  about  355°.  These  seismograms  have  signals  with  much  higher  frequency 
than  those  shown  in  the  other  two  seismograms  recorded  at  an  azimuth  of  about  313°.  Thus,  Q 
around  Garm  station  is  perhaps  more  complicated.  We  expect  to  establish  this  in  our  next  phase 
of  waveform  modeling. 

Comparison  of  Reeional  PnJ  Waves  from  the  US  and  USSR  Crustal  Models:  In  this  section,  we 
continue  to  investigate  the  regional  waveforms  that  are  likely  to  be  predicted  by  the  crustal  models 
developed  for  North  America  and  Soviet  Union.  Since  the  record  at  Harvard  station  was  so 
successfully  modelled  and  since  it  was  at  a  range  of  640  km,  we  examined  the  response  of  the 
USSR  crust  model  at  this  range.  In  fact,  seismograms  at  such  distances  are  just  becoming  available 
from  the  Soviet  Union.  Figure  46  shows  the  comparison  between  the  synthetic  seismograms 
computed  for  the  samples  of  the  USSR  and  US  crustal  models  for  both  the  vertical  and  radial 
components  assuming  the  same  focal  mechanism  at  the  same  azimuth.  The  seismograms  were 
computed  at  the  respective  depths  of  the  earthquakes  and  the  two  depths  are  similar.  We  also 
used  the  same  source  function.  The  GARM  crustal  model  predicts  a  stronger  PmP  relative  to  the 
Pp.  The  Sn/SP  and  sSn  arrivals  predicted  in  the  seismograms  by  the  US  crustal  model  (marked  by 
the  arrows)  seem  to  exhibit  a  correspondence  to  the  long-period  signals  predicted  in  the  response 
of  the  USSR  crustal  model.  The  difference  in  the  amplitude  ratios  is  caused  by  the  differences 
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Comparison  between  broad-band  vertical  component  data  at  GARM  station  and 
A I  ,  5*'^  seismograms  computed  for  different  Q-models  including  the  elastic  model. 

Also  shown  ,s  the  first  75  seconds  of  each  seismogram  for  each  Q  model. 
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Figure  45.  Broad-band  seismograms  recorded  at  GARM  station  from  several  azimuths  where  the 
variation  of  frequency  content  observed  on  these  seismograms  is  a  clear  effect  of  crustal  Q  along 
the  propagation  paths. 


Comparison  of  Oisplocement  Seismogroms  ot  640  Km 
for  the  U.S.  and  USSR  Crustal  Models 
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Figure  46.  Comparison  between  the  synthetic  seismograms  computed  at  640  km  for  the  USSR 
crustal  model  and  the  US  crustal  model.  Both  vertical  and  radial  components  are  shown. 
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in  the  near-surface  velocities  of  the  two  crustal  models.  For  the  North  American  crustal  model, 
the  crustal  velocities  have  a  gradient  near  the  surface.  The  rays  arrive  at  the  receiver  more  steeply 
compared  to  the  rays  for  the  USSR  crustal  model,  thus  partitioning  the  energy  in  a  significantly 
different  ratio  to  the  vertical  and  radial  component. 

We  further  investigated  the  composition  of  the  Pnj  and  Sni  waves  predicted  by  the  USSR 
crustal  model  at  640  km  in  terms  of  generalized  rays.  We  found  that  the  S-wave  reflections  from 
the  Moho  and  Conrad  discontinuities  are  strong  as  shown  in  the  top  two  seismograms  of  Figure 
47.  The  Conrad  reflection,  ScS,  arrives  immediately  following  the  Moho  reflection  SmS.  The 
phase  shown  by  an  arrow  on  the  SmS  seismogram  is  arriving  at  the  arrival  time  of  S„  phase,  but 
its  waveshape  is  more  complicated  than  is  expected  from  a  classical  Sn  phase.  Among  the  other 
phases  that  contribute  most  significantly  to  the  total  P„|  seismogram  within  the  S  window  are  the 
sSmS,  SmSSmS  and  ScSScS  phases.  The  bottom  two  seismograms  plotted  in  Figure  40  allow  us  to 
compare  the  agreement  between  the  generalized  ray  (Total)  and  frequency  wavenumber  (F-K) 
seismograms.  The  agreement  is  poor  following  the  PmP  arrival.  A  strong  long-period  signal  does 
propagate  to  the  receiver  in  the  frequency  wavenumber  seismogram.  This  must  be  a  total  effect 
of  many  generalized  rays.  This  effect  was  also  observed  on  the  recorded  seismogram  at  Garm 
station  even  at  a  distance  of  200  km  from  the  source. 

Conclusions:  Based  on  these  investigations,  it  seems  feasible  to  develop  time-domain  discriminants 
at  different  nuclear  test  sites  which  rely  on  the  stable  features  that  are  observed  in  the  recorded 
waveforms.  The  most  stable  phases  are  observed  in  the  explosion  generated  Ppi  waveforms  for 
periods  as  short  as  2  s  (Burdick  ei  al.,  1991c;  Saikia  and  Burdick,  1991).  In  this  study,  we  have 
extended  our  analysis  approach  to  regional  broadband  seismograms  from  earthquake  sources.  The 
short-period  P^i  waves  have  a  functional  dependence  on  the  crustal  waveguide.  They  can  be 
deterministically  modeled  using  average  flat-layered  crustal  structures  and  using  some  selected 
generalized  rays.  By  modeling  the  broadband  displacement  at  Harvard  station,  we  found  that  the 
structure  across  the  crust-mantle  transition  zone  and  within  the  mantle  can  profoundly  affect  the 
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Figure  47.  Understanding  the  waveform  computed  at  640  km  from  the  USSR  crustal  model  using 
the  ray  decomposition  technique.  The  top  six  seismograms  are  for  the  individual  ray  groups. 
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frequency  content  of  the  phases  like  S„  and  sS„.  The  source  multiplicity  of  an  earthquake  can  also 
create  added  complexity  in  the  frequency  content  of  these  phases.  We  found  that  the  seismograms 
near  the  S-wave  arrival  can  adequately  be  modeled  using  the  ray  responses  of  the  following  phases: 
SmS,  sSmS,  SmSSmS,  sSmSSmS  and  SmS'SmS. 

For  the  Soviet  Union,  the  most  important  requirement  for  understanding  recorded  seismo¬ 
grams  is  the  crustal  model.  The  structure  within  the  Soviet  Union  is  heterogeneous  and  the 
development  of  reliable  crustal  models  is  on-going  (Gurrola  and  Minster,  1991).  In  this  study,  we 
have  developed  a  crustal  structure  by  modeling  the  recorded  seismograms  at  Garm  station  from  an 
earthquake  at  a  distance  of 200  km  ( Az=292°)  which  consisted  of  Conrad  and  Moho  discontinuities. 
In  addition,  a  slight  gradient  is  allowed  for  the  upper-mantle  structure.  The  ray  analysis  indicated 
that  the  most  important  generalized  rays  for  the  composition  of  the  s„i  waves  are  the  following 
phases:  SmS,  ScS,  sSmS,  sScS,  SmSSmS  and  sSmSSmS. 

We  found  that  the  phases  with  longitudinal  propagation  mode  is  less  affected  by  varying 
attenuation  structure.  The  phases  like  sSmS,  SmSSmS  which  travel  mostly  with  shear-mode  of 
propagation  are,  however,  significantly  affected.  Therefore,  the  spectra  of  the  shear-wave  phases 
will  be  deleted  in  high-frequency  and  the  spectra  ratio  of  PmP  and  PmPPmP  phases  to  sSmS, 
SmSSmS  phases  may  provide  a  better  event  discriminantion.  For  the  problems  related  to  the 
estimation  of  yield  of  a  nuclear  explosion,  the  dependence  of  the  waveform  amplitude  on  the  Q 
models  can  be  crucial.  It  can,  however,  be  seen  in  Figure  38  that  the  intial  Pnl  waveform  (say,  first 
20s  after  the  Pn  onset)  is  less  sensitive  compared  to  the  Snl  waveforms  to  the  chosen  variation  in 
Q  models.  Thus,  it  may  appear  that  the  analysis  based  on  the  early  part  of  the  regional  Pnl  may  be 
a  better  source  for  estimating  yield  of  nuclear  explosions.  The  explosions  are  expected  to  be  shallow. 
If  they  are  buried  within  the  low  Q  material,  it  is  likely  that  the  effect  may  be  severe.  Further 
investigations  are  needed  to  find  these  effects. 


88 


0  STUDIES  PART  II  -  PRELIMINARY  RESULTS  FROM  BROADBAND  MODELING 
OF  LONG  RANGE  REGIONAL  SEISMOGRAMS 
Introduction:  Recent  studies  in  the  United  States  suggest  that  a  great  deal  of  information  about 
structure  and  attenuation  be  obtained  from  braodband  modeling  of  regional  phases  as  demonstrated 
in  the  previous  section  (Saikia  and  Burdick,  1992).  In  this  pilot  effort,  we  have  performed  some 
basic  calculations  to  test  the  usefulness  of  present  earth  models  in  predicting  the  observations. 
Figure  48  displays  a  map  of  Asia  and  digital  stations  recording  two  events  from  the  Hindu-Kush 
region.  Since  many  events  of  different  focal  mechanisms  occur  in  this  region  at  various  depths,  it 
becomes  an  ideal  source  region.  Two  events  have  been  studied  in  detail  as  reported  by  Zhao  and 
Helmberger  (1991,  see  Figure  49).  Broadband  seismograms  at  several  upper  mantle  distances  are 
modeled  in  these  studies  using  an  earth  model  appropriate  for  the  Tibet  region,  see  Figure  50.  The 
three  stations,  namely  K1 V,  OBN  and  ARU,  lie  towards  the  northwest  of  this  source  region  and  the 
earth  structure  to  the.se  stations  appears  to  be  the  most  homogenous.  In  Figure  51  through  52,  we 
show  broadband  recordings  of  two  events  of  July  24, 1989  (3h  27m  48.77s)  and  February  5, 1990 
(5h  16m  45.0s)  at  these  sites.  Both  the  events  have  simple  seismograms  with  similar  looking  S  at 
ARU.  The  P-waves  are  the  most  dissimilar,  having  different  polarities  and  different  strengths  for 
pP  and  sP  depth  phases.  Both  events  produced  motions  that  rotate  well  into  (P-S V)  and  (SH).  The 
same  is  true  for  the  other  stations  although  the  noise  level  at  OBN  is  panicularly  strong.  The  P-waves 
are  nearly  nodal  for  the  1990  events  at  these  stations. 

From  the  previous  studies  of  shield  regions  in  North  America,  we  would  expect  the  paths  to 
ARU,  OBN  and  KIV  to  be  similar  to  the  paths  decribed  by  SNA  model  for  S  waves  (Grand  and 
Helmberger,  1985)  and  by  S25  model  for  P  (LeFevre  and  Helmberger,  1989),  see  Figure  50. 
However,  we  have  a  problem  in  the  source  region  where  the  crustal  thickness  is  much  greater  than 
in  SNA,  roughly  60  km  vs.  35  km.  TJius,  the  broadband  synthetics  predicted  by  these  models  (SNA 
&  S25)  show  an  obvious  problem  in  that  sP  and  sS  occur  too  early,  see  Figures  51  and  52.  Note 
that  observed  sS  is  especially  late  on  the  1990  KIV  record.  Since  this  event  is  at  a  depth  of  102 
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Figure  49a.  Teleseismic  modeling  of  broadband  P- waves  assuming  a  simple  surface  interaction, 
namely  pP  and  sP,  which  fixes  the  source  depth  and  allows  relocation  (Zhao  et  al.,  1991).  (Event 
7/24/89) 
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Figure  49b.  Teleseisiuic  modeling  of  broadbsnd  P->vaves  assuming  a  simple  surface  interaction, 
namely  pP  and  sP,  which  fixes  the  source  depth  and  allows  relocation  (Zhao  et  al.,  1991).  (Event 
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Figure  51.  Vertical  Component:  broadband  comparison  of  synthetics  and  observations  where  the 
synthetics  were  generated  with  reflectivity,  restricted  to  periods  greater  than  1  sec. 


94 


Figure  52.  Radial  Component;  broadband  comparison  of  synthetics  and  observations  where  the 
synthetics  were  generated  with  reflectivity,  restricted  to  periods  greater  than  1  sec. 
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kms,  we  would  expect  this  problem  to  occur.  Another  computational  difficulty,  at  present,  is  the 
inability  to  include  attenuation  directly  into  our  synthetics.  The  intrinsic  Q  structure  for  a  given 
crustal  medium  can  be  handled  in  the  frequency-wavenumber  code.  This  code  has  been  calibrated 
as  discussed  in  the  previous  section  and  will  be  used  in  future  studies  to  monitor  the  effect  of  Q  on 
broadband  signals.  Alternatively,  we  can  simulate  the  appropriate  behavior  by  convolving  these 
results  with  a  t*  operator.  For  a  world-wide  average  value,  a  t*  of  1  for  P  waves  and  a  t*  of  4  for  S 
waves  are  normally  used.  The  later  procedure  is  adopted  in  this  present  investigation. 

A  comparison  of  these  synthedcs  with  observations  indicate  some  agreements  and  some 
disagreements.  For  example,  the  synthetic  S  waves  show  good  agreement  with  the  recorded  data 
at  most  stations.  The  P  wave  synthetics  show  some  inconsistencies  with  the  data,  especially  at 
ARU  89  where  the  polarity  is  even  wrong.  But  this  much  disagreement  is,  in  overall,  expected 
given  that  the  source,  structure  and  t*  are  not  known  well  at  this  stage  of  modeling. 

An  enlarged  portion  of  the  SV-wave  for  the  best  looking  station,  ARU,  is  given  in  Figure  53 
along  with  various  t*  operators.  A  t*  between  2  and  3  appears  to  fit  the  waveshape  the  best  although 
some  adjustments  in  the  triplications  are  needed  to  improve  the  level  of  fits. 

It  is  generally  useful  to  decompose  the  synthetics  into  rays  so  that  individual  pulses  can  be 
isolated  and  modeled,  and  the  agreement  between  the  data  and  the  synthetic  is  improved.  Figures 
54  and  55  display  the  synthetics  containing  the  simple  surface  reflections,  pP,  sP,  and  pS,  sS.  These 
three  arrivals  do  quite  well  at  matching  the  reflectivity  synthetics  at  the  nearest  distances  but  less 
well  as  the  distance  increases.  This  is  cau.sed  by  the  neglect  of  the  S-P  interactions  in  the  crust  or 
the  so-called  S  coupled  PL  waves.  A  detailed  comparison  of  the  S-waves  at  OBN  as  generated  by 
the  reflectivity  code  indicates  that  a  small  P-wave  precursor  occurs,  see  vertical  and  radial.  This 
is  caused  by  the  SV-to-P  conversion  at  the  crustal  receiver  structure.  This  feature  has  not  been 
included  in  the  rays  at  this  stage.  The  precursor  in  the  ray  synthetics  are  caused  by  diffracted 
P-waves  along  the  crustal-mantle  interface. 
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Figure  53.  Comparisons  of  observations  with  reflectivity  synthetics  at  ARU  indicating  the  promise 
of  detailed  fits  with  some  adjustments  in  upper- mantle  triplications. 
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Figure  54.  Vertical  Component  comparison  of  synthetics  with  observations  when  the  phases  P,  pP, 
sP  and  S,  pS,  sS  have  been  generated  by  ray  theory  assuming  a  t*  =  1  and  4. 
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Figure  55.  Radial  Component  comparison  of  synthetics  with  observations  when  the  phases  P,  pP, 
sP  and  S.  pS,  sS  have  been  generated  by  ray  theory  assuming  a  t*  =  1  and  4, 


Included  in  the  ray  synthetics  is  the  absolute  amplitudes.  In  many  cases  the  synthetics 
assuming  t*=I  and  4  bracket  the  data,  at  least  for  the  event  with  the  best  mechanism,  namely  the  90 
event.  The  timing  of  the  surface  reflected  phases  in  these  synthetics  is  different  than  is  the  reflectivity 
run  because  we  used  deeper  sources  to  compensate  for  the  thicker  crust,  that  is  h=102  and  1 17  kms 
respectively. 

It  is  difficult  to  check  the  timing  of  P  and  S  in  these  figures  because  of  the  mismatch  in 
amplitudes.  These  features  are  more  easily  seen  in  synthetics  where  the  radiation  patterns  have 
been  supressed  as  in  Figures  56  and  57.  These  figures  contain  only  the  down-going  P  and  S  where 
the  sources  are  at  their  proper  depths,  85  and  106  km  respectively.  Only  the  vertical  componant  is 
displayed  which  yields  the  clearest  arrival  onsets.  We  have  included  pure  tectonic  stlye  synthetics 
for  comparison.  The  P-waves  are  slightly  early  for  the  S25-SNA  model  but  the  SV  waves  are  early 
by  10  secs  or  more,  see  ARU  and  KIV.  On  the  other  hand,  from  the  GCA-TNA  model  the  SV 
waves  are  late  at  KIV  (2  secs)  and  ARU  (8  secs).  The  travel  times  are  in  general  agreement  with 
the  (SS-S)  times  reported  on  by  Woodward  and  Masters  (1991 )  and  Grand  and  Helmberger  (1985), 
see  Figure  58.  Results  from  the  latter  study  indicate  a  sharp  increase  in  velocity  when  crossing  the 
high  ridge  of  topography  that  extends  from  Hindu-Kush  to  Lake  Baikal.  Velocities  north  of  this 
feature  appear  unusually  fast,  essentially  SNA.  This  study  as  well  as  the  earlier  report  by  Rial  et 
al.  (1984)  used  the  waveform  and  travel  times  of  upper  mantle  (SS-S)  data.  The  results  from 
Woodward  and  Masters  ( 1 99 1 )  are  teleseismic  type  (SS-S)  measurements  and  can  contain  anomalies 
deeper  than  Upper-Mantle.  Whatever  the  reason  they  see  much  more  variability  in  velocity 
structures  across  Russia  than  previously  concluded,  along  the  Ural  mountain  extension  from  NZ  to 
the  Hindu-Kush.  A  deep  anomaly  running  North  and  South  beneath  the  Mid  United  States  has  also 
been  seen  by  Grand  (1991),  perhaps  this  is  a  similar  feature. 
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Figure  58.  SS-S  residuals  in  Eurasia.  Negative  residuals  are  indicative  of  faster  than  average 
velocity  material,  while  positive  residuals  indicate  slower  velocity  material,  after  Woodward  and 
Masters  (1991) 
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The  most  interesting  events  occur  along  the  tectonic  regions  indicated  by  the  topography 
display  and  are  probably  in  slower  structure  than  most  of  the  Soviet  stations.  Thus,  to  model  the 
many  broadband  records  available  would  require  approximateing  these  structures  by  at  least  2  D 
models. 

A  common  method  of  generating  syntheucs  for  2D  models  is  due  to  Chapman  (1978),  called 
the  WKBJ  method.  Some  short-period  results  for  models  formed  by  a  linear  connection  between 
pure  models  is  given  in  Figure  59,  see  Helmberger et  al.  ( 1 985).  Five  profiles  are  displayed  showing 
the  upper-mantle  triplications.  The  small  first  arrival  at  ranges  near  1 8°  is  essentially  a  lid  diffraction 
coming  from  a  depth  of  about  185  kms  in  the  pure  S25  model.  The  large  second  arrival  at  18°  for 
this  model  is  just  the  "400"  triplication.  Since  the  lid  velocities  are  slower  in  the  TIP  model  a  third 
arrival  coming  from  a  depth  near  250  km  is  apparent  int  he  pure-TlP  synthetics.  The  other  three 
profiles  correspond  to  placing  the  source  at  various  positions  along  the  1000  km  transition  zone. 
Profile  S25-T1P  is  appropriate  for  a  midpoint  and  the  other  two  are  mostly  SN  A(-500)  or  TIP(+500). 
The  synthetics  for  (S25-TIP  fits  the  absolute  travel  times  at  ARU  quite  well.  The  code  used  in 
constructing  these  synthetics  have  been  used  primarily  in  modeling  multibounce  SH  waves  (Grand 
and  Helmberger,  1985;  Graves  and  Helmberger,  1988)  and  discussed  in  detail  in  previous  WCC 
reports  (WCCP-R-83-01). 

A  new  technique  of  generating  synthetics  for  2D  structures  is  presently  going  through  the 
testing  process  but  will  allow  broadband  modeling  since  diffractions  and  tunneling  will  be  treated. 
This  method  is  essentially  an  extension  of  the  Cagneard  de-Hoop  method  to  handle  lateral  variations. 
Figure  60  displays  the  basic  geometry  and  rayset  for  a  direct  P  or  S.  Following  this  approach 
requires  finding  a  geometrical  ray  parameter  (po)  or  snell’s  law  angle  (0o)  that  tracks  rays  f’‘om  the 
source  to  the  recevier  after  reflecting  from  each  boundary.  The  shallowest  generalized  rays  con¬ 
tribute  both  headwaves  and  post  critical  angle  reflections.  The  timing  of  these  various  rays  can 
also  be  used  to  generate  synthetics  directly  as  suggested  by  Burdick  and  Salvado  (1986). 
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Figure  59.  Profiles  of  WKBJ  synthetics  (P-waves)  assuming  2D  structure.  Panels  S25  and  TIP  are 
pure  path  results  where  the  P-waves  are  about  3  secs  faster  for  SNA  relative  to  TIP.  The  synthetics 
for  a  model  that  starts  with  TIP  and  ends  with  S25  fits  the  timing  at  ARU  quite  well. 


TL 


S25 


,<&S^4?^»i««W///////' 


m*»00MHt  tft- 
*mmm0m»0ati0  iirri^i^ 
m^^00000»f*  •nr  '^itt 

—  900m000f0‘tlfi>rrt^ 


km\\\\\\\^^NNC)X«»s^ 
l.'m\\\\\\\\VN«NKNHOB« 
~  mW^^SWKKIi^AIKMN!* 
.m\\\\\\\\\\\VMimiMi 
Km\\\\\\\\\\KM»M!03 


'^^MW/zz/mmiJiirir 

K0tMM0/t/t/mnmf 

^0!MM/*m/mtrttt 
rjK0im/0//i//iiwiinr 


vmmsiti«s<siaseiaieiBK^00'00^^ 
~mm\\\\\\\\\\»iM«ieK!a^^ 
^vmmsiKiixta»siKK0W0M/fmHmir 
-vmsssfi^s99SBOoeiaKm!000fwmmM 
-immiiS90ssssaear^j!M0f//Mmmr 

^^nseoseaaaarMgMf/iM/mtm 
vmiissskBeeee00Mnrfr/rnmm~ 
- womr . . . 


-vssmiiii!isaes0S0M0000m/m/r 

mmxiiiS9sesisw00M//mim 

■Kv\\\\\v\\v^mo!wr/Ka////M 

mk\\\\\\\\K\xnr000z//y//mm 

in^\\\\\N\v>mim/////^H 

mmk\\\\\\\^jr0/’d////£mm 

mmmiiWMxmsw/mmmmm 


kvxwa^Amm/M 
-k\xwi\r/snxMr~ 
miKscmmam 
mewiwjstxam 
— kiaw/xms 


7  11  0  10  n  0 

Vrlocit  V 


-100  GOO 

Distnncc  (  ktn  ) 


n  9  10  1 

vclocil  V 


gure  60.  Display  of  ray  paths  connecting  a  source  and  receiver  to  reflectors  used  in  generating 
nthetics  for  2D  models. 


In  conclusion,  it  appears  that  an  excellent  estimate  of  regionalized  attenuation  can  be  con¬ 
structed  from  broadband  data.  However,  to  obtain  high  resolution  will  require  refining  the  pro- 
pagational  corrections.  This  appears  possible  by  applying  the  latest  analytical  techniques  to  the 
numerous  data  sets,  namely  explosions  (Gamero  et  al.,  1992)  and  earthquakes  (Zhao  et  al.,  personal 
communication). 

Conclusions  and  Recommendations:  The  release  of  internal  data  from  the  Soviet  Union  obviously 
presents  some  great  opportunities  for  solving  some  of  the  most  long-standing  yield  estimation 
problems.  However,  interpretation  of  the  data  is  still  complex,  and  we  continue  to  face  the  challenges 
of  modeling  the  source  RDP,  determining  its  dependence  on  depth  and  separating  source  from 
propagation  effects.  Here  we  have  focused  on  the  influence  of  Q.  We  recommend  that  integrated 
Q,  velocity  and  source  models  be  developed  for  the  Soviet  Union  which  yield  consistent  values  for 
regional  and  teleseismic  t*.  The  Soviet  PNE  program  has  been  much  more  extensive  than  the  U.S., 
and  intensive  studies  of  this  data  should  yield  important  new  insights  into  source  scaling.  Wc  also 
recommend  an  intensive  effort  to  obtain  near  field  explosion  data  from  the  Soviet  Union  since  data 
of  this  type  has  proved  so  valuable  in  interpreting  U.S.  data. 
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